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ABSTR4CT 


A numerical investigation of two-dimensional natural convection flow 
and heat transfer from a substrate-mounted flush heat source immersed in 
a liquid-filled square enclosure was conducted. The study is relevant to 
direct liquid-immersion cooling of electronic components. A control 
volume based finite-difference model that accounts for conduction heat 
transfer within the substrate and heat source and the coupled natural 
convection in the fltiid was utilized. Numerical predictions were obtained 
for a wide range of Rayleigh and Prandtl numbers, substrate to fluid and 
heat source to fluid thermal conductivity ratios and other geometrical 
parameters that may be encountered in practice. An increase in Rayleigh 
number lead to more vigorous flow and promoted cooling. No noticeable 
effect on the nondimensional temperatures was observed when changing 
the Prandtl number from 7 to 100. Little reduction in maximum 
temperatures was observed when substrate and component to fluid thermal 
conductivity ratios were increased beyond 10 and 25, respectively. 
Component to substrate width ratio change from .25 to .999 resulted in 
approximately linear decrease in the maximum temperature. A 
companion experimental study of three-dimensional natural convection 
transport from a flush mounted array of heat sources in water was also 
conducted. Computed temperatures compared favorably to appropriate 
experimental data. 


m 





TABLE OF CONTENTS 


4 


I. INTRODUCTION .1 

A. STATEMENT OF PROBLEM.1 

B. OBJECTIVES.4 

II. NUMERICAL MODEL . 7 

A. MATHEMATICAL FORMULATION. 7 

B. METHOD OF SOLUTION.9 

III. PARAMETRIC STUDY .12 

A. EFFECT OF Ra.12 

B. EFFECT OF Pr.24 

C. EFFECT OF SUBSTRATE CONDUCTIVITY RATIO R,.24 

D. EFFECT OF COMPONENT CONDUCTIVITY RATIO Rc.41 

E. EFFECT OF COMPONENT TO SUBSTRATE WIDTH 

RATIO w/dt.51 

IV. EXPERIMENT .61 

A. EXPERIMENTAL APPARATUS.61 

1. Additions.61 

2. Changes.62 

B. LIQUID CRYSTAL.63 

V. NUMERICAL COMPARISON WITH EXPERIMENT.73 


A. WHAT TO COMPARE ?.73 

B. 0.2 WATT NUMERICAL COMPARISON.73 

C. HIGHER POWER NUMERICAL COMPARISON.80 


IV 

























D. UPPER/LOWER HEATER NUMERICAL EVALUATION.82 


VI. CONCLUSIONS.91 

VII. RECOMMENDATIONS.93 

APPENDIX AUSER PORTION OF FINITE DIFFERENCE PROGRAMS...94 

APPENDIX B CONTOUR PROGRAM.106 

APPENDIX C TRAVERSE PROGRAM.116 

APPENDIX D TEMPERATURE ACQUISITION PROGRAM.117 

APPENDIX E UNIFORM POWER PROGRAM.119 

APPENDIX F UNIFORM DATA REDUCTION PROGRAM.122 

APPENDIX G NONUNIFORM POWER PROGRAM.125 

APPENDIX H NONUNIFORM DATA REDUCTION PROGRAM.128 

LIST OF REFERENCES.132 

INITIAL DISTRIBUTION LIST.134 



















NOMENCLATURE 


Cp specific heat at constant pressure [J/kg-KI 

d dimension in Fig. 1 [m] 

dt substrate thickness [m] 

g gravitational acceleration [m/s^ 

h component height [m] 

H enclosure height [m] 

k thermal conductivity [W/m-K] 

dimension in Fig. 1 [m] 

p pressure [N/m^] 

P non-dimensional pressure, p/pU^^ 

Pr Prandtl number, pcj/kf 

q non-dimensional heat fiux [eq.7] 

Q Heat generation rate per unit length [W/m] 

P Power to the heater 

Rc protrusion to fluid thermal conductivity ratio, k^/kf 
R, substrate to fluid thermal conductivity ratio, k^ 

Ra Rayleigh number, gPQhS/ak^v 

51 non-dimensional counter clockwise contour distance along the 

component-fluid/component-substrate interface in Table 2 

52 non-dimensional coxmter clockwise contour distance along the 

solid-fluid interface in Table 1 


VI 




Sb, St dimensions in Fig. 1 [m] 

Lead resistance [ A ] 

Rp Precision resistor resistance [ A ] 

% enclosure wall temperature [K] 

T non-dimensional temperature, (t-t^VCO/kf) 

u vertical velocity component [m/s] 

U non-dimensional vertical velocity component, u/Uq 

Uq reference velocity, (gPQh/kf)i^ [m/s] 

V horizontal velocity component [m/s] 

V non-dimensional horizontal velocity component, v/Uq 

Vl Voltage measured over the leads and heater [Volts] 

Vi< Voltage over the precision resistor, leads and heater [Volts] 
w component width [m] 

X vertical coordinate [m] 

X non-dimensional vertical coordinate, x/h 

y horizontal coordinate [m] 

V non-dimensional horizontal coordinate, y/h 

Greek Symbcls 

a fluid thermal diflusivity [m^/s] 

P coefficient of thermal expansion [1/K] 

5X non-dimensional vertical distance between grid point and the 
control volume face 

5Y non-dimensional horizontal di8tan<» between grid point and the 
control volume face 


vu 



p 

p 


fluid density [kg/in^] 
dynamic viscosity [kg/m-s] 

V kin ematic viscosity [m^/s] 

K nondimensional conductivity [R/(RaPr)i^] 

Skibacripts 

c component (chip) 

f fluid/liquid 

max maximum 

s substrate 

sur substrate or component surface 
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I. INTRODUCTION 


A. STAT 


Wi’j I w 


NT OF PROBLEM 


Heat removal from electronic equipment is an area of extensive research 
and will be more important in the future due to the ever increasing 
component volumetric heat generation rates. The advantages of liquid 
cooling have been demonstrated in the thermal control of macroelectronic 
components such as power transistors and power supplies for almost 40 
years. However liquid cooling of microelectronic equipment has only 
recently gained wide attention due to the need for increased power 
dissipation, the availability of inert dielectric liquids, and the introduction 
of immersion cooled supercomputers. Bergles and Bar-Cohen [Ref. 1] and 
Nakayama [Ref. 2] have described various cooling techniques currently 
under investigation to meet modem component cooling requirements. 
Incropera [Ref. 3] discusses convective cooling. 

In natural convection the flow is induced by buoyancy forces. The 
buoyancy force is generated when a body force acts on a fluid in which a 
density gradient is present. Natural convection cooling has the advantages 
of simplicity of design, low operating cost due to no outside power 
requirements and minimum maintenance, absense of noise and high 
reliability. Also, in the event of a mechanical failure, natural convection 
may be the only method of cooling available [Ref. 4-5]. Natural convection in 
liquids results in much higher cooling rates compared to those of air and 
hence numerous studies have recently been performed on natural 
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convection from discrete heat sources inside liquid filled enclosures [Ref. 6- 

103 . 

Gaiser [Ref. 6] experimentally investigated natural convection liquid 
cooling of a three column array of uniformly heat generating simulated 
electronic components on a flat plate and in a vertical channel by using foil 
heaters. The channel was made by placing a smooth movable shroud 
parallel to the high density heater surface. Temperature measurements 
and flow visualizations at various power levels and channel spacings were 
made. Experimental correlations relating local heat transfer with shroud 
wall spacing and foil heater power levels were developed. Flow was 
laminar for all power levels studied. For the unshrouded case, the flow of 
coolant over the center column resulted in a weak fluid entrainment from 
the side columns and in turn enhanced the cooling of the side columns. 
For the shrouded case, the fluid entrainment by the center column was 
much greater and had a greater cooling effect on the side columns. For a 
channel of 3 mm or less, the increased frictional resistance caused a 
significant increase in the component temperature and the correlations 
were no longer followed. 

Joshi, Willson and Hazard [Ref. 7] conducted an experimental 
investigation of steady state and transient natural convection from a 
column of rectangular heated protrusions in a vertical channel in water. 
Flow visualizations and component surface temperature measurements 
were carrried out for several power dissipation levels and channel 
spacings. For the smallest channel spacing, the component surface 
temperatures increased significantly due to a reduction in fluid velocity. 
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During transient periods, an initial diffusive transport was observed 
followed by the evolution of convective effects. No overshoots in component 
temperatures were observed. Steady transport responses over several 
component power levels were also examined. 

Park and Bergles [Ref. 8] simulated microelectronic circuits with flush 
mounted and protruding thin foil heaters. Heat transfer coefficients were 
obtained for various heater heights and widths in water and R-113. The 
heat transfer coefficient was found to increase with decreasing width, 
which was attributed to three dimensional boundary layer effects. 
Protruding heaters were found to have a coefficient about 15 percent above 
flush mounted heaters. Coefficients for upper flush heaters were found to 
be less than lower heaters. For protruding heaters, the upper heaters had 
higher coefficients than the lower heaters. 

Bar-Cohen and Schweitzer [Ref. 9] provided an analytical basis for the 
design and optimization of convective immersion cooling systems. 
Analytical development and experimental verification of the relations for 
the natural convection heat transfer coefficients prevailing along the 
surfaces of immersed, uniformly heated plates in both symmetric and 
asymmetric configurations were conducted. Milanez and Bergles [Ref. 10] 
experimentally studied the wake temperature distribution produced by a 
horizontal line source on a vertical adiabatic surface and compared it to the 
similarity solution. They also examined the interaction of two heaters. 
Velocities induced by the lower heater enhanced the heat transfer 
coefficient at the upper heater, however this effect was offset by the 
preheating of the fluid by the lower heater. Jaluria [Ref. 11] conducted a 
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numerical study of heat sources mounted on an adiabatic wall under the 
influence of natural convection. Sathe and Joshi [Ref. 5] conducted a 
numerical investigation of natural convection flow and heat transfer 
arising from a substrate mounted protruding heat source. Their study 
accounts for conduction heat transfer within the protrusion and substrate 
and the coupled natural convection in the fluid. 

Areas of current concern with the new packaging techniques are stress 
on the electronic components due to thermal expansion coefficient 
mismatch. Temperatures throughout the package are needed so thermal 
stress analysis can be conducted and high stress areas can be located. The 
complex geometries of the boards and the coupled conduction/convection 
heat transfer problem lead to intricate heat transfer paths. The difficulty is 
magnified when the circuit board spacing is small and the flow pattern is 
populated by numerous non-uniform heat sources. An examination of the 
heat transfer and fluid flow characteristics of such a system would aid 
designers in the development of better cooling schemes. 

B. OBJECTIVES 

This study was performed to determine the effects of changing the 
various parameters involved in the construction of electronic packages on 
the heat transfer and to observe the effects of a noniiniform high density 
grouping of discrete heat sources in a vertical channel for direct liquid 
natural convection cooling. Gaiser [Ref. 6] describes the high density 
grouping of heat sources. 
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The study consisted of two parts. First a computer model was utilized to 
conduct a two dimensional study to find the effects of various governing 
parameters on the natural convective heat transfer and fluid flow behaviors 
for a substrate mounted flush component immersed in a liquid. Then the 
three dimensional flow over a simulated printed wire circuit board 
consisting of multiple arrays of flush-mounted heaters was observed 
experimentally. Results of the computer model were compared to results 
from past and current experimental studies. 

Specific objectives of the numerical parametric study were: 

- To determine temperatures and flows throughout and around a 
simulated direct liqtiid natural convection cooled flush moimted electronic 
component. 

- To study the effect of various combinations of parameters for an 
electronic component flush mounted on a substrate. 

- To investigate the effect of two component interaction on the heat 
transfer and flow. 

The experimental study had the following objectives: 

- To determine the effect of channel spacing on fluid flow and heat 
transfer. 


5 





- To observe the efTect of nonuniform powering on heat transfer and fluid 
flow. 


• To observe the surface temperature patterns during uniform and 
nonuniform heating using liquid crystal. 
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11. NUMERICAL MODEL 


A. MATHEMATICAL FORMULATION 

The schematic diagram of the configuration examined is shown in Fig. 
1. A flush heat source mounted in a vertical substrate is immersed in a 
two-dimensional fluid-filled square enclosure of height H. The dark shaded 
region represents the heat source and the light shaded area the substrate. 
Uniform volumetric heat generation takes place in the flush mounted 
heater. The enclosure boundaries are maintained at a constant 
temperature, t^. The heater, the substrate and the fluid have constant but 

different thermophysical properties. Assuming a steady-state, laminar 
flow with no viscous dissipation and the Boussinesq approximation to be 
true, the dimensionless governing equations can be written as: 


FUtidRegian 


9U ^ 

ax dY 


= 0 


( 1 ) 


a(u2) aojv) 
ax av 


ax2 av^ ox 


(2) 


d(uv) ■ acv^) 
ax av 




a(UT). a(VD 
ax av 




(3) 

(4) 
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Figure 1. Schematic diagram ofoonfigoration for numerical study 













SoUd Region (heater) 

h 

ax2 av^' wRc (5) 


SoRdBegion (eubetrate) 


a^ ^ a^ 
ax2 av^ 


= 0 


( 6 ) 


where X=x/h, Y=y/h. U=uAJo. V*v/Uo, T=(t-te)/(Q/kf), P=p/pUo2; 
Uo*(gpQh/kf)i^, Pr*|iCp/kf, Rasg^Qh^/ok^v, R^sk^/kf and R,sk^. 

The boundary conditions constitute isothermal enclosure walls at 
temperature t^ and the no-slip and impermeable wall conditions for the 

velocity components; i.e UsVsTsO at the enclosure walls. Heat fluxes are 
appropriately matched at the interfaces of dissimilar materials. The 
following governing parameters emerge as a result of the non- 
dimensionalized governing equations, boundary and matching conditions: 
Ra, Pr, Rc, R,, Sj/h, d/h, Lt/h, S^/h and w/d^ The definitions of the symbols 

can be found in the Nomenclature. 


R METHOD OF SOLUTION 

The governing equations are discretized using a finite difference scheme 
wherein the control volumes for the temperature and pressure are 
staggered from those for the velocities. Power law profiles are used for the 
spatial variation of the dependent variables to ensure realistic results for a 



wide range of the grid Peclet numbers. Interface diffusivities are 
calculated using a harmonic-mean formulation in order to handle abrupt 
changes in the material properties. The details of the discretization process 
can be found in Patankar [Ref. 12]. 

The discretized equations were solved iteratively using the line-by-line 
TDMA (Tri-Diagonal Matrix Algorithm) and the SIMPLER procedure as 
outlined in Ref. 12. It is noted that even though separate equations are 
written for the fluid and solid regions, the solid region is ntunerically 
simulated by letting its viscosity become very large. Thus the same 
momentum equation is solved throughout the computational domain. 
Similarly, for the energy equation, heat sources and property values are 
used implicitly in the control volume formulation wherein mass, 
momentum and energy balances are incorporated for individual control 
volumes [Ref. 12]. 

Test computations were performed on a series of grids ranging from 
30x30 to 73x73 control volumes to determine the grid size effects. The values 
and locations for the maximum temperature and peak velocities did not 
change appreciably when the grid was refined beyond 40x40 control 
volumes. Most calculations were performed on a grid size of 50x50. 
However for Ra > 10^ the solution would not converge on a 50x50 grid and it 
was found that by increasing the control volumes in the area of the 
suspected high thermal gradient, convergence could be achieved. 
Therefore for Ra s 10^ the number of control volumes was increased to 
70x70. Convergence was based on a balance of the rate of energy generated 
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in the component and the rate of energy leaving the enclosure walls. 
Convergence was based on the energy balance being better than 1 percent. 

The program was run on the Mechanical Engineering VAX station 2000 
Cluster and the IBM-370/3033 computer. The greater the number of control 
volumes the the greater the time to run on each system. Also if Ra was of 
the order of 10^, it took up to 120,000 iterations for covergence to take place. 
The relaxation factor for the velocities, temperatures, pressure and buoyant 
force were adjusted until a proper combination was found that would make 
the solution converge. After coming up with a proper set of relaxation 
factors, they could be used for other high Ra solutions and the number of 
iterations reqtiired was reduced. The run times could also be reduced 
when previously computed solutions were used as an initial guess. The 
VAX system took approximately 10 times longer to nm a program with the 
same number of iterations and control volumes than the IBM system. 
However, more VAX stations could be used on a continuous basis which 
made using the VAX system slightly preferable. 





m. PARAMETEUC STUDY 

Numerical computations were conducted for various Ra, Pr, R,, and R^ 
values that would be found in actual cooling applications [Ref. 5]. Ra 
numbers from 10^ to 10^ were investigated with the limit of 10^ due to the 
very slow convergence of the computations for Ra greater than 10^. The 
values of Pr, R, and R^ were varied from 7 to 100, 0.1 to 100 and 1 to 100 

respectively. The geometric parameters used in the computations were 
Sb/hs2, d/hs:2, Lt/h=2, St/h=3, and w/dt varied from 0.25 to 1. A baseline 
case with Ra=10^, Prs25, R^slO, Rcs25 and w/d^sO.S was conducted and 
when a particular parameter was being studied all other parameters were 
kept the same as in the baseline case. Table 1 shows the breakdown of the 
cases studied. Appendix A shows the user portion of the finite difference 
program used to obtain the numerical data. Appendix B contains the 
contour program used in making the following figures. 

A. EFFECT OF Ra 

The temperature contours for Ra^lO^to 10^ are shown in Figs. 2(a)-(d). 
These are nearly symmetric with respect to the enclosure vertical 
centerline for all Ra. The high R, (10) value tends to spread the heat out to 

each side even though the heater is placed on the right side of the substrate. 
As Ra increases, the extent of thermal stratification increases in the lower 
portion of the enclosure. All boiindary layers get thinner and a larger 
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Table 1. iwirameters (aP panuneters except the one being studied 

revert to these) and the range of parametCTS 


Parameter 

Baseline 

value 

Range studied 

Ra 

105 

103, 104^ 105^ 106 

Pr 

25 

7,25,100 

R. 

10 

.1,1,10,50,100 

Rc 

25 

1,25,50,100 

w/dt 

.5 

.25, .5, .75, .999 


thermal gradient occurs in the fluid next to the substrate and the enclosure 
top and sides. 

Figures 2(a)-(d) also show changes in the streamline patterns with more 
vigorous flow with an increase of Ra. The flow travels up on both sides of 
the package. The flow forms two cells with the right cell being somewhat 
larger and stronger, due to the direct heater-fluid interaction. Flows are 
nearly S 3 rmmetric on the right and left and as Ra is increased the eddy 
centers move up and toward one another. Also as Ra is increased, the fluid 
in the lower right portion of the enclosure becomes increasingly stagnant. 

The heat loss though each face of the component is shown in Fig. 3. As 
Ra is increased the percentage of the heat lost directly to the fluid i.e. from 
the front face, is increased due to the increased convection. As Ra 
increases the percentage of heat lost from the top and bottom of the 
enclosure decreases as shown in Fig. 4, and the heat lost to the sides of the 
endosure increases. The heat loss through the endosure bottom is below 5 
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Figure 4. Heat loss throu^ midosure walls for various Ra 













percent for all Ra shown but it decreases with an increase in Ra due to 
increased thermal stratification near the bottom. There is a greater 
percentage of heat loss through the left wall than the right even though the 
right wall faces the component and the flow is greater on the right. This is 
because the left wall is closer to the substrate and heater than the right 
wall. It is interesting to note that the top wall dissipates almost half the 
heat generated by the component. 

Figure 5 shows the substrate nondimensional surface heat flux (q) 
which is defined as: 



Table 2 describes the values of S2 corresponding to its location on the 
substrate. Heat fluxes are not defined at the comers. S2s0 at the lower 
right comer of the substrate and goes counter-clockwise around the solid- 
fluid boimdary including the component-fluid boundary. The same general 
trend for the variation of heat flux with S2 is seen for all Ra. The heat flux 
goes through a minumum at S2»0.5 after which it rapidly reaches its 
maximum value at S2=2, at the base of the component. The heat flux then 
decreases going up the component and when the substrate above the 
component is reached at S2s3 it decreases even more steeply, going through 
a minimum at S2<«4.5 prior to reaching the top right comer. At the top 
right corner the heat flux decreases sharply again going over the top 
substrate face and comes to the minimum heat flux on the substrate at the 





Table2. Ranges of S2o(nTesp(nidii]£ to dififerentsolid-fluidinterfiaoes 


Surface no. 
in Fig. 1 

X 

Y 

S2 

1 

0to2 

0 

0to2 

2 

2to3 

0 

2to3 

3 

3to5 

0 

3to5 

4 

5 

Oto-1 

5to6 

5 

5to0 

-1 

etoll 

6 

0 

-ItoO 

11 to 12 


midpoint of the top face at S2s5.5. Reviewing the temperature contours in 
Figs. 2(a)-(d), this is where the fluid temperature is the hottest next to the 
substrate. The heat flux on the top face is nearly symmetric. Proceeding 
down the substrate left face the heat flux goes through a small dip but 
quickly increases as one approaches the back of the component. It reaches 
a local maximum heat flux at the back of the base of the component at S2 « 9 
and then decreases to a local minimum at S2sl0.25. It then increases and 
comes to a local maximum at the bottom left comer at S2=ll. The heat flux 
on the bottom face shows the same symmetry as the top, however it is much 
higher due to colder fluid in the vidnity of the bottom face. 

Nondimensional substrate surface temperatures (T,ur) for various Ra 

are shown in Fig. 6. The substrate surface temperature is computed using 
the harmonic mean formulation as; 


T = 

^ fttr 




( 8 ) 
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The lowest temperatures are at the bottom of the substrate and the highest 
are at the top of the component. Local maxima are evident at the back of the 
substrate due to direct conduction from the heater. The nondimensional 
temperatures are higher for the lower Ra. The higher the Ra, the greater 
the buoyancy induced flow and hence higher heat transfer from the 
component, leading to lower non-dimensional temperatures. Note that if Q 
and kf are held fixed, a higher T implies a higher actual (dimensional) 

temperature. In such a case, higher Ra can be obtained by using a liquid 
with higher p, leading to higher buoyant forces, enhanced convection and 
reduced actual temperatures. Thus the dimensional results are consistent 
with the non-dimensional trends. 

Figure 7 shows the nondimensional component surface heat flux for 
various Ra versus SI. Table 3 shows the corresponding position of the 
component to SI with SlsO at the lower right comer of the component and 
increasing going cotmter-clockwise around the component. Fluxes for all 
Ra show similar trends. The component front face, in direct contact with 
the fluid, has the lowest heat flux. The top face has the second highest heat 
flux and the back face has a slightly lower heat flux. The bottom of the 
component has generally the greatest heat flux. The largest heat flux is 
near the bottom right comer. 

The trends of Fig. 8 show the component nondimensional surface 
temperatures to be similar for all Ra. Surface temperatures are fairly 
xiniform arotmd the component with the highest temperature near the top 
on the front face. The high Rg and R« values allow the heat to be uniformly 

spread so that the component temperature remains fairly uniform. 
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Maximum temperature versus Ra is shown in Fig. 9(a). A decrease in the 
maximum temperature from 0.122 to 0.042 is observed for an increase in Ra 
from 10^-10^. Fig. 9(b) shows that the change in maximum temperature 
versus Ra is similar for all w/d^. Although nondimensional temperature 

decreases, the dimensional temperature may go up because the non- 
dimensional temperature is inversely proportional to the heat transfer 
coefficient. An increase in Q would bring about an increase in the 
dimensional temperature but a decrease in the nondimensional 
temperature. As Q increases, the natural convection flow velocities 
increase and this would increase the heat transfer coefficient. 

REFFECTOFPt 

Fluorinert dielectric liquids used in immersion cooling have Pr values 
much greater than 1. Computations for Pr=s7, Prs26 and PrslOO were 
therefore conducted. These revealed that with other parameters remaining 
constant, Pr had a negligible effect on temperature contours, streamlines, 
heat flux or temperatures as is seen in Figs. 10-17. This would indicate that 
the results of experimental studies in water with Prs7 would apply to the 
dielectric liquids as well for the same value of Ra and other parameters. 

C. EFTECT OF SUBSTRAIE CONDUCnvnY RATIO Rg 

The streamlines and temperature contours for various R, values when 
RaslO^ are shown in Figs. 18(a)-(e). For .1^R,^1 all the main flow is 

generated by the heater side with a very weak cell on the left side of the 
enclosure. When R, reaches 10 both cells are of about the same strength 

and the streamlines form two cells that start on opposite sides of the 
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Figure 10(a). Temperature oontoiira and streamlines for Prs7 
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Figure 10(b). Temperature contours and streamlines for PrB25 
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Figure 11. Heat loss through oomponent faces for various Pr 
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Figure 12. Heat loss throu^ endosiune walls for various Pr 
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Figure 18(b). Temperature contours and streamlines for R,«l 

















































Figure 18(d). Temperature oontoun and streamlines for R^sSO 
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substrate. Temperature contours and streamlines for 10 100 

essentially look the same and it appears that there is little change in the 
induced flow beyond R,=10. Temperature contours show a very high 
gradient just outside the component for R^sO.! and the hot spot of the 
component is in the top left comer. For R^sl the hot spot moves towards the 
center of the component and as R, increases the component hot spot moves 
towards the fluid. 

Figure 19 shows the percent heat loss through the component faces for 
various R. with Ra=10S. There is a tremendous decrease in the percent 
heat loss to the fluid fmm the front face when R, is increased from 1 to 10 
after which it shows only a slight decrease. As the conductivity of the 
substrate increases it is easier for the heat to flow through faces other than 
the front. The substrate acts like a fin, increasing the surface area where 
heat loss can occur. The back face percent heat loss increases throughout 
the range of R, increase. The bottom and top face percent heat loss increase 

when Rs is increased from 0.1 to 10 and then remain fairly constant. 

The percent heat loss through the enclosure walls for various R, is 

shown in Fig. 20. The bottom wall initially shows a slight decrease when 
R, increases from 0.1 to 1 due the left cell of the flow in the fluid bringing 

higher temperatures closer to the bottom wall and then a slight increase 
over the remaining range of R, due to the increased heat conduction 

through the bottom portion of the substrate into the vicinity of the bottom 
wall. The other percent heat loss changes are a combination of the 
increased thermal conductivity and flow changes. 
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Figure 21 shows the substrate surface heat flux for various R,. For 
R,^10 the heat flux for the substrate is fairly uniform. For R,<10 most of 
the heat flux comes from the front face of the component for S2 ranging 
from 2-3. The greatest heat flux is at the base of the component at S2s2. 
The same pattern follows for the substrate surface temperatures in Fig. 22 
with the temperatures being uniform for R,^10 and for R,<10 the 

temperatures being higher over the component surface, S2s2-3, with the 
highest temperature at the top of the component, at S2s3. 

Component surface heat flux is shown in Fig. 23. For R,<10 the heat 

flux through the component-fluid interface is greatest, for Sl^O to 1. For 
Rg^lO the component-substrate interface shows the higher heat flux with 

the bottom face having the maximum. Very little difference occurs when 
R, is increased beyond 50. The component surface temperature is fairly 
uniform for all R, due to the high conductivity of the component as shown in 
Fig. 24. The surface temperature decreases greatly for R,s0.1 to 10 and 

then decreases slightly over the remaining range with very little change 
between R,=50 to 100. As the substrate conductivity increases there is less 

resistance to the flow of heat and hence the surface temperatures decrease. 

Maximum temperature versus R, is shown in Fig. 25 and shows there is 
little benefit to increasing R, beyond 10. 

D. EFFECT OF COMPONENT CONDUCnVITY RATIO R* 

The effects of R^ on the temperature contours and streamlines are 
shown in Fig. 26 (a)-(d). For R^sl there is a steep thermal gradient in the 
component with the maximimi temperatiire near the center of the 
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Figure 22. Substrate surface tenq w rat u re for various R^Ra«10» 
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Figure 26(d). Temperature contours and streamlines for RckIOO 
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component. For Res25 the hot spot is towards the fluid face and the 
component is almost at a uniform temperature. For R^sSO and 100 the 

temperature contours look identical. There is a negligible change in the 
streamlines for the complete range of R^ studied. 

Figure 27 shows the percent heat loss through the component faces. For 
Rc increasing from 1 to 25 there is a large decrease in the percent heat loss 

from the front face and a slight decrease from the back. A large increase in 
the percent heat loss from the bottom and a slight increase in the top face 
make up for the decrease at the front and back. With the increase in 
component conductivity, heat flows easier through the component and a 
greater percent can then go out the top and bottom. The substrate has a 
relatively high conductivity and the top and bottom subst ite areas act like 
fins. For R^ increasing from 25 to 100 there is a slight increase in the 

bottom percent heat loss and a very slight decrease from the front and the 
back. The bottom face uses the substrate below as a fin to move the heat to 
the lower cooler region of the enclosure. The heat loss through the 
enclosure walls, shown in Fig. 28, changes insignificantly which is 
expected due to no change in the fluid flow or the temperature distribution 
in the fluid over the entire range of R^ studied. 

The component surface heat flux. Fig. 29, shows for R^sl a similar 
pattern for each face. Each face has a local maximum flux at the center 
with the largest maximum flux at the center of the back face at Sls2. The 
low conductivity of the component limits the movement of heat, so the 
greater flux is associated with the surfaces on the component that have the 
most volume of component aroimd them. These locations on the component 
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Figure 28. Heat loss tfaroui^ emdosure walls for various R^, Ra^KK^ 
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surfaces correspond to the center of each face. For 25 100, the 

component surface heat fluxes are nearly identical with the greatest flux 
going through the bottom face near the fluid-component boundary. The 
flux is greatest here since the substrate extending below the component 
bottom surface acts like a fin. In addition there is induced flow over the 
substrate below the bottom of the component augments heat transfer from 
the front substrate surface below the heater. 

The component surface temperatures for R^sl vary widely with the 
maximum occuring at SI *0.7 as seen in Fig. 30. For 25 ^R^^ 100 the 

temperature over the surface is fairly uniform. The maximum 
temperature versus R^ in Fig. 31 shows little advantage in increasing R^ 

higher than 25. 

E. EFFECT OF COMTONENT TO SlJBSrRATEWlDlH RATIO w/dt 

Streamlines and the temperature contours in the liquid region in Figs. 
32(a)-(d) show only minor changes when w/dt changed from 0.25 to 0.999. 
This occurs because the high R, and R^ result in low solid thermal 
resistance, with the major temperature drop occuring in the fluid. In the 
component region the hot spot moves from the fluid surface to the 
component center as w/dt increases from 0.25 to 0.999. 

The percent heat loss through the component faces is shown in Fig. 33. 
For w/dt increasing from 0.25 to 0.999 there is approximately a 10% 

increase of the heat loss from the top face and a 20% increase from the 
bottom face. Both the top and the bottom percent heat losses rise because of 
the increased heat transfer area as w/dt increases and a greater volume of 
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Figure 32(a). Temperature oontours and streamlines for w/dt«0i!5 
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Figure32(b). Temperatureoontonrianditreamlinwfor 
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Figure 32(d). Temperature contours and gtreamlineB f6r w/dt«a999 
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the component next to them. There is approximately a 7% decrease in heat 
loss {torn the front face because more heat is being taken away from the top 
and bottom and the front is actually becoming less of the total surface area 
as w/dt is increased. There is approximately a 23% decrease on the back 

face. The back face now is also less of the total percent area of the 
component surface and the back has less of the substrate acting as a fin as 
w/dt increases until there is essentially no substrate at w/d^ = 0.999. There 

is no difference in the percent heat loss from the enclosure boundaries, Fig. 
34, as was suggested from the temperature contours and streamlines. 
Figure 35 shows the maximum temperature versus w/d^. There is 

almost a linear decrease in maximum temperature with an increase in 
w/d^. The same rate of heat is being dissipated by a larger volume which 

reduces the maximum temperature. 


58 







PERCEENT HEAT LOSS 



■ FRONT 
0 TOP 
B BACK 
E BOTTOM 


Figure 33. Heat loss through component faces for various w/d^ RaslQs 
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IV. EXPERIMENTS 


A. EXP] 


D4 NI 104 


^AL APPARATUS 


1. Additions 

The details of the experimental apparatus can be found in Gaiser 
[Ref. 6]. It consists of a test surface with 3 columns of flush heaters with 15 
heaters per column as shown in Fig. 46. The test surface is immersed in a 
large plate glass tank filled with deionized water. Each heater has a 
thermocouple attached to the midpoint of the back. Voltages of the 
thermocouple and heater leads are taken by a data aquisition system and a 
HPCHewlett Packard) computer determines the temperatures and power of 
the heaters. Further additions and changes done in the experimental setup 
are described below. 

A 3-dimensional traverse assembly with controller for the traverse 
motors and a deaei .ting system for the water were installed. The traverse 
controller interfaces with the HP computer in the two horizontal directions 
and a program was developed to move the traverse with the HP computer 
integrating the temperature acquisitions with a probe. The traverse must 
be manually moved in the vertical direction. See Appendix C for the 
traverse computer program. Deaeration was achieved by spraying water 
through a nozzle into a partially evacuated bottle. The water was then 
allowed to accumulate under the vacuum and pumped back into the tank at 
a rate similar to which it was being drawn out. Teflon sheets were placed 
over the tank's free surface to reduce the aeration of the water after it had 
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gone through this process. When the deaeration system was activated and 
a slight (3" or greater Hg) vacuum was drawn, the bottle would collapse 
leaving the system unoperable. A nonflexible container would be desirable 
for a better operation of the deaeration system. 

2. Changes 

A change was made in the way the power input to the heaters was 
calculated. Upon taking over the apparatus it was noticed that the 
resistances of the various heaters varied by as much as 8% from the lowest 
resistance reading to the highest. The high density heater board was taken 
apart and the heaters with the resistances farthest from the norm were 
replaced. The new heaters were of a different batch than the original 
heaters and they all had a resistance of 10.4 ohms originally but when they 
were measured through the leads after they were installed the resistances 
again varied by as much as 6%, with the lowest value being the same as the 
heater before installation and the highest resistance being 6% above the 
original heater resistance. Since all the heaters were the same before being 
installed it was assumed that the excess resistance came from the leads. 
The excess resistance for all leads was then determined by subtracting the 
original heater resistance from the resistance of the heater and the leads 
and a computer program that took this into account was implemented. 1% 
and 5% precision resistors in parrallel with the original 1% precision 
resistor were used to increase the cizrrent of the heaters that had less power 
due to the excess lead resistance. The following relation was used to 
determine the correct power: 
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V J (9) 

where the symbols are identified in the Nomenclature. Each heater and 
lead resistance had to be treated individually and in some cases extra line 
resistance had to be added to the leads to make the heaters have all the 
same power. A uniform and nonuniform power program were developed 
for the HP computer. The uniform power program was used when each 
heater was given the same power. The nonuniform power program was 
used when the power to the side columns was different than the center 
colimm of heaters. Also nonuniform and uniform data reduction progams 
were developed to obtain the temperatures and to reduce and store the data 
simultaneously for each case. See Appendixes D-H for computer programs 
utilized in the acquisition of experimental data. 

a UQUID CRYSTALS 

Hallcrest thermochromic liquid crystals whose optical properties of 
the planar texture give rise to the selective reflection of bright colors that 
change with temperature [Ref. 13] were used in the experiment. The 
liquid crystals had to be kept refngerated prior to usage. The first step in 
liquid crystal application was to obtain a dean dry surface. A black backing 
was then applied with an air brush. After the black backing had dried, the 
liquid crystals were mixed with a binder with 3 parts binder to 1 part liquid 
crystal slurry and the mixture was applied over the backing using the air 
brush. Care was taken to keep the layers very thin to reduce any 
interference with the heat transfer process. The water in the tank had to be 
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kept very clean and the water filtering system was run for two days 
constantly before placing the test surface with the liqmd crystals on it in the 
water. If any contamination was present, only a dark blue color would 
show up which was difficult to see on the black backing. Even with the 
water being very pure, the liquid crystals deteriorated and could only be 
expected to last less than two weeks. It is believed that the crystals 
somehow react with the polar water molecule and become inactive. In 
dielectric fluid this had not been a problem, which was evident in past 
thesis research [Ref. 14]. 

The crystals were calibrated against temperatures using a platinum 
resistance standard (Rosemont Engineering Model 920A Commuting 
Bridge) and a Model 913A Calibration Bath. Two ranges of liquid crystal 
sensitivity to temperature were used. After standardization for the colors 
which appeared below 30 degrees C, the crystals compared within 0.5 
degrees C with the thermocouple temperatures of the heaters, but at higher 
temperatures (greater than 40) it was difficult to determine the change 
from blue to black and an agreement between the crystals and the heater 
thermocouples was about 5 degrees C with the thermocouple temperature 
always higher than the color indicated. 

A problem with working with the crystals was to photograph the 
resulting surface temperature color patterns. It is very difficult to shoot 
into a glass container filled with water. All light sources outside the tank 
must be turned off or all that will be seen in the picture is the refection of 
the photographer. The correct light source was also an important 
consideration. Florescent lights do not have a complete spectrum of light 
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being emitted and the pictures taken under these conditions did not show 
the colors that were actually present. Taking pictures with a spot lamp on 
the face with a high speed film (ASA 1600) did not work because the area 
with the spot ended up being burned out in the picture and the surroimding 
area of the picture appeared black. The best method found was to obtain a 
light source that contained the wavelength band necessary. It was 
obtainable at a local photograph store and it was used with a film speed of 
400 ASA or higher. Care had to be taken to keep the light source above the 
water and shining inside the tank. Shielding was placed so the bulb would 
only light up inside the tank staying away from the refecting glass and 
lighting only the test surface. The light was only switched on when taking 
pictures to avoid setting up fluid variations by warming up the tank water. 

Figure 36 shows the view that is shown in the pictures to follow. It also 
shows the shroud and the test surface. The shroud is a dear plastic that 
causes a constriction to the flow and simulates another circuit board or a 
wall in an actual situation in electronic equipment. At no power, the test 
surface appears black and corresponds to ambient temperature (18 degrees 
C). The first picture, Fig. 37 is a picture with no shroud and a power level 
of .5 watts per heater. The actual colors that were present are brown which 
appears at 24.0 degrees C, green which appears at 25.3 degrees C, and blue 
which appears at 28.3 degrees C and dissappears at 44.0 degrees C where 
the surface returns to black. The black and white pictures can only be used 
to show higher temperature regions with the lighter regions on the 
photographs indicating where the hot spots are. The color patterns become 
more evident near the middle and top compared to the bottom of the heater 
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Figure36. Schematic of the view that is seen in the IbBowing pictures 



66 













columns. The temperature patterns seem symmetric for all three columns 
and there appears to be no thermal interaction between the three columns. 
Heat being brought up by natural convection is brought up the face and 
some is transfered back into the plexiglass substrate. As the power is 
increased with no shroud the color patterns become more evident at the 
bottom of the heater column. At three watts per heater power input, Fig. 38 
shows again that there is no interaction between the columns. The color 
pattern above the top heater indicate higher temperattires in that region 
compared to the 0.5 watts per heater case. 

A shroud is placed 6mm in front of the test surface in Fig. 39 and 
entrainment of the flow is evident by the shape of the temperature patterns. 
The shroud prevents entrainment flow from the front and water must come 
from the sides towards the center column. This actually enhances cooling 
of the outside top heaters for a shroud spacing of 6 mm as seen in Figs. 39 & 
40, compared to when there is no shroud. The bottom heaters appear not to 
be effected by the shroud placement. In Fig. 41 the shroud spacing is 3 mm 
and the power input per heater is 0.2 watts. The color patterns now become 
visible at a much lower power level. The flow at the top is greatly 
constrained with the high temperature patterns actually merging over the 
center column. Temperatures are increased drastically at the 3mm shroud 
spacing and Fig. 42 shows that at a 0.5 watt power input per heater the 
temperatures are as high as with a 1 watt, 6 mm shroud spacing case or a 2 
watt, no shroud case. 

Nonuniform power levels for the side columns and center column at 
different power input were looked at for no shroud and 3mm shroud cases. 


67 












Figure 38. No shroud case, 3 watts per heater power input 



Figure39. 6 nun shroud spacing,^ watt per heater power input 
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Figure 40. 6 mm shroud spacing, 1 watt per heater power input 



Figure 41.3 nun shroud spacing, J2 watt per heater power input 
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Figs. 43-45 , and showed the same patterns as the uniform power cases. 
For the no shroud case there did not appear to be any interaction between 
the three columns. For the small shroud spacing cases there was 
considerable side entrainment. It was difficult to go to higher power levels 
in the low shroud spacing cases because air bubbles formed over the top 
center region and disrupted the flow and made the temperatures increase 
greatly even for power levels of the order of 0.5 watts per heater. 





ii 


Figure42. 3 mm shroud spadng, US watt per heater power input 
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Figure 44. No shroud case, 0.5 & 1 watt per heater power input 

(oenter/sides) 
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Figure 45. 6 mm shroud spacing, .2&.5 watt per heater input (center/side) 
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V. NUMERICAL COMPARISON WITH EXPERIMENT 


A WHAT TO COMPARE ? 

The liquid crystal temperature patterns show that there are pronoimced 
3-dimensional efifects in the flow and temperature patterns with low shroud 
spacings, especially near the top of the 3 columns of heaters. With no 
shroud and wider shroud spacings, the heat transfer and fluid flow 
problem is primarily 2-dimensiona?, especially at the lower center portion 
of the 3 columns and at the center of the heaters. The central plane of the 
bottom heater of the center coltimn (section AA in Fig. 46) was chosen for a 
2-dimensional numerical comparison. 

R 0^ WATT NUMERICAL COMPARISON 

The first attempt to model the temperature measured with a 
thermocouple at the middle of the back of the heater for the 0.2 watt case 
agreed only by approximately 25%. Figure 47 shows the heater surface 
temperatures versus SI. SI, as in the parametric study, is the counter¬ 
clockwise contour distance aroimd the heater starting at the lower right 
comer of the heater. The experimental value of the thermocouple is located 
at SI slightly greater than 1.5. To get closer to the experimental 
temperature the numerical model was examined again. A two heater 
model was developed to account for the effect of upper heaters. Figure 48 
shows the corresponding temperature contours and streamlines. With the 
two heater model the agreement was within approximately 13%, as seen in 
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SECTION AA 



Figure 46. Schematic diagram of the experimental heater studied 
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Figure 47. experimental/numerical compi 
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Fig. 49. The heaters were placed within the grooves in the plexiglass test- 
surface using a low thermal conductivity epoxy. During the liquid crystal 
work it was observed that the epoxy thermal conductivity was lower than 
that of plexiglass, so, a thermal conductivity of y2 that of plexiglass was 
used in the numerical model to see what effect it would have on the 
comparison. It proved to be a slight improvement over the plexiglass 
thermal conductivity case but not as much of an improvement as the two 
heater model, as seen in Fig. 50. The dimensions of the heater were then 
looked at closely. The heater is made up of a Kapton sheet and an inconel 
sheet together. The dimensions of the Kapton were the dimensions of the 
heater that were being used but the inconel was actually the resistor that 
was causing the increase in temperature and its dimensions were smaller 
than the Kapton. The height of the region which contained the inconel was 
0.0062 m and the width was 0.022 m. When these new smaller dimensions 
were used in the numerical model the agreement was much better as seen 
in Fig. 51. Temperatures for both thermal conductivities (0.14 W/m-K and 
0.07 W/m-K) are plotted in Fig. 51 and the agreement with the experiment is 
within approximately 9% for the one heater model. 

The heater was a layer of Kapton over a sheet of inconel as stated above. 
When the heater was modeled as two layers, the numerical prediction was 
within 7 % of the experimental measurement as shown in Fig. 52. When 
the two heater model was used with tiie smaller dimensions for the heater 
and the heater divided into two layers, as mentioned above. Fig 53 shows 
that the temperature variation for the two different substrate conductivities 
bracketed the experimental result, as was expected. 
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Figure 49. Pint attempt at 2 heater modd 















Figure 52. One heater nuxtel, heater split into a Kapton and Inconel layer 











C. HIGHER POWER NUMERICAL COMPARISONS 

The measured component temperatures at higher power levels 
experimental values did not agree as closely as the 0.2 watt power case. 
Figure 54 shows that as the power input per heater is increased the 
disagreement between the numerical and experimental values increases. 
At 0.5 watt power input the percent difference is approximately 8 %. One 
possible reason for this is that as the power is increased, the downstream 
heaters have a greater effect on the upstream heater temperatures due to 
greater conduction through the substrate. Thus the experimentally 
measured temperatures are higher by a greater extent compared to the 
numerical model. Thus for comparison with higher power cases, 
additional heaters must be added both on the side and on the top in the 
model. Another possible reason for the percent difference increase with 
power is that the properties for the water used in the model were taken at 
the average temperature of the heater and the ambient water. Also, the 
inconel region of the heater was taken as a uniform heat source in the 
model but in reality is a ribbon that snakes up and down on the Kapton 
sheet. At higher power levels this could cause hotter areas near the ribbon 
and this could cause the numerical value to be slightly lower. Still another 
possibility is that there are 3-dimensional effects in reality that the'2- 
dimensional model cannot detect and these effects get larger with increased 
power. Entrainment of the warmer fluid from the sides columns into the 
center region of the heater covdd make the experimental value greater than 
the 2-dimensional model. 
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EXPERthCNTAL VALUE 
SIBSTRATE K..14 W/m-K 
SUBSTRATE K«.07 W/m-K 




D. UPPER/LOWER HEATER NUMERICAL EVALUATION 

The numerical model for two heaters provided the upper and lower 
heater temperatures. A study of the effects of the lower heater on the upper 
heater was conducted. Figures 55-59 show the progression of the top and 
bottom heater temperatures as the power is increased from 0.2 to 3 watts 
equally in both heaters. Figures 60-61 collect the temperature patterns at 
the various power levels studied for the two heaters. Similar patterns are 
observed as power level increases for all top heaters and all bottom heaters 
when compared with themselves. When comparing the lower heater to 
upper heater temperatures, there is a larger increase in the top heater's top 
face temperatures. The back face temperatures do not decrease as much on 
the top heater as the bottom heater. Also, the bottom face temperatures of 
the top heater do not decrease below the temperature of the bottom of the 
front face contrar:^ to the bottom heater. It is more difficult for the upper 
heater to lose heat because of the natural convection bringing up the heat 
from the lower heater, thus the upper heater temperatures are higher than 
the lower heater temperatures. The front face temperatures of both heaters 
are always closer than the other faces. As power level increases, the top, 
back and bottom faces become farther from their counterparts on the bottom 
heater. 

Figures 62-65 compare the heat loss through the heater face for 
increasing power levels. The only face that shows an increase for the top 
and the bottom heater is the front face. At low power levels the percent heat 
loss through the front face with increasing power is its largest. Both 
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Figure 55. 0.2 watt case top and bottom heater surfiEU» temperatures 



Figure 56. 0 l 5 watt sur&oe temperatures for the top and bottmn heater 
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Figure 63.0.2^ watt top face percent heat loss for top and bottom heater 



Figure 64.0.2 *3 watt back &oe percent heat loss for top and bottom heater 
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heaters show a similar pattern. The back and bottom faces percent heat 
loss both decrease for increasing power levels 

The effect on the lower heater of a heater above it is shown in Figs. 66-67. 
For 0.2 watt power input per heater shown in Fig. 66, the one heater model 
surface temperatures are all less than the corresponding 2 heater model 
siuface temperatures for the bottom heater For 3 watts power per heater, 
shown in Fig. 67 the opposite is true. The surface temperatures for the 
bottom heater of the 2 heater model are lower except at the top of the heaters 
where it is equal to the one heater model. The increased convection for the 2 
heater model when the heat input is much higher causes lower 
temperatures at the bottom of the lower heater. Thus for 0.2 watts, 
conduction from the upper heater to the lower heater causes an increase in 
temperature of the bottom heater, whereas for a higher power input of 3 
watts, the increased convection caused a cooling effect on the bottom heater. 








Figure 65.0^ watt bottom face percent heat loss for top and bottom heater 



Figure66.0.2 watt surfooe temperatures fcnrbottcnn heater for 1&2 heater 

models 
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1 HEATER MODEL 

2 HEATER MODEL 










VL CONCLUSIONS 


During the course of the numerical and experimental study, several 
important conclusions were drawn about the transfer of heat fium a flush 
heat source. A constant flux of temperature boundary condition must be 
used carefully as the numerical study showed that it can be far from 
uniform when using the conjugate solution to solve the problem. With the 
advent of faster computers it seems that all heat transfer problems should 
be solved with a combined conduction/convection procedure to improve 
accuracy of the solution. 

Changes in the Pr between 7 and 100 did not change the nondimensional 
temperatures. This was extremely important to this study because the 
experimental work has been done in water, with a Pr of 7 and most fluids 
that are currently being used for the inunersion of electronic components 
are fluorinerts and have a Pr greater than 10 and less than 100. The 
experimental results should therefore be applicable to these fluids as well. 

Increasing the R, provides substantial cooling enhancement when Rg is 

less them 10 but increasing it past 10 only provides slight enhancement. For 
the designer of electronic equipment, this should help in the decision as to 
how high a conductivity the substrate should be. When a choice of more 
than one material is available, the cheapest material with an R, of 10 or 

greater could be chosen. 

An increase in R^ provides lower and more uniform component 
temperatures primarily for R^ less than 25. The designer can once again 
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look for the cheapest material that would do the job with an of 25 or 
greater and optimize the cooling of the component as well. 

There is almost a linear decrease in maximum temperature with 
increase in w/d^. This would suggest that the componet should be made as 

wide as the substrate if the design permits. 

Finally, liquid crystals proved to be a valuable tool in determining the 
areas of the hot flow up the front of an experimental surface. The effect of 
the shroud spacing became very evident as the spacing decreased. 
Although at small shroud spacing the 3-D effects of the fluid flow were 
evident, at large spacing and no shroud conditions and near the bottom of 
the high density board always the heat transfer process could be modeled as 
2-dimensional. Very good agreement of numerical and experimental 
results were seen, especially at the low power levels. 
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VII. RECOMMENDATIONS 


In continuation of this study, it is suggested that the same configuration 
of three high density heater columns be used, however it is recommend that 
the following aspects be fiirther investigated: 

-Utilize liquid crystals to determine the spacing where shroud effects are 
negligible 

-Use a 3 heater model to determine what effects this would have on the 
lower heater at the higher power levels 

-Analyse the effects of using non-steady power inputs on the heat 
transfer 

-Utilize a 3-dimensional numerical model to see if 3-D effects are 
significsmt, especially at the high power levels 

-A nonflexible container should be acquired for the dearation system to 
allow the use of the test surface at higher power levels 

-High power studies to determine the onset of turbulence could be looked 
at with the accompanying effects on the heat transfer 
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APPENDIX A 

USER PORTION OF FINITE DIFFERENCE PROGRAMS 


c **********tHESIS********************•’> *************************** 

c 

BLOCK DATA 

LOGICAL LSOLVE,LPRINT,LBLK,LSTOP 

COMMON F(75,75,5),P(75,75),RHO(75,75),GAM(75,75),CON(75,75), 

1 AIP(75,75),A1M(75,75),AJP(75,75),AJM{75,75),AP(75,75), 

2 X(75),XU(75),XD1F(75),XCV(75),XCVS(75), 

3 V(75),yv(75),YDIF(75),YCV(75),YCVS{75), 

4 YCVR(75),YCVRS(75),ARX(75),ARXJ{75),ARXJP(75), 

5 R(75),RMN(75),SX(75),SXMN(75),XCVI(75),XCVIP(75) 

COMMON DU(75,75),DV(75,75),FV{75),FVP(75), 

1 FX(75),FXM(75),Fy(75),FYM(75),PT(75),QT(75) 

COMMON/INDX/NF,NFMAX,NP,NRHO,NGAM,Ll,L2,L3,Ml,M2,M3, 

IIST,JST,ITER,LAST,TITLE(13),RELAX{13),TIME,DT,XL,YL, 

2IPREF,JPREF,LSOLVE(10),LPRINT(13),LBLK(11),MODE,NTIMES(10), 
3RHOCON,ZERO 
COMMON/CNTL/LSTOP 
COMMON/SORC/SMAX,SSUM 
COMMON/COEF/FLOW,DIFF,ACOF,ALFA 
DIMENSION U(75,75),V(75,75),PC(75,75) 

EQUIVALENCE(F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3), 
1PC(1,1)) 

DIMENSION TH(75),THU(75),THDIF(75),THCU(75),THCVS(75) 
EQUIVALENCE(X(1),TH(1)),(XU(1),THU(1)),(XDIF(1),THDIF(1)), 

1(XCV(1),THCV),(XCVS(1),THCVS(1)),(XL,THL) 

DATA NFMAX,NP,NRH0,NGAM/5,6,7,8/ 

DATA LSTOP,LSOLVE,LPRINT/1*.FALSE.,10*.FALSE.,13*.FALSE./ 

DATA MODE,TIME,ITER/1,0.,0/ 

DATA RELAX,NTIMES/13*!.,10*1/ 

DATA LBLK/11*.TRUE./ 

DATA DT,IPREF,JPREF,RHOCON/l.E+10,1,1,1.0/ 

C ENTER DATA INTO LABELED COMMONS USING DATA STATEMENTS 

C****«******************************************«****************** 

DATA TITLE(1) ,TITLE(2),TITLE(3),TITLE( 4) , 

1TITLE(6)/'VELU','VELV','STFN','TEMP' , 

2'PRES'/ 

DATA (LSOLVE(I),I-l,4),LSOLVE(6),LPRINT(3)/6*.TRUE./ 

DATA RELAX(1),RELAX(2),RELAX(4),RELAX(6),ALFA/.04,.04,.4,.4,.!/ 
DATA LAST/5000/ 

C DATA LAST/5/ 

END 

(^*************************************************«********** 
SUBROUTINE USER 

^**************************************«****«**************** 


94 









COMMON F(75,76,5),P(75,75),RH0(75,75),GAM(75,75),CON(75,75), 

1 AIP(75,75),AIM(75,75),AJP{75,75),AJM(75,75),AP(75,75), 

2 X(75),XU(75),XDir(75),XCV(75),XCVS( 75) , 

3 Y(75),yv(75),YDIF(75),VCV{75),yCVS(75), 

4 YCVR(75),YCVRS(75).ARX(75),ARXJ(75),ARXJP(75), 

5 R(75),RMN(75),SX(75) ,SXMN(75),XCV1(75),XCVI P ( 75) 

COMMON DU(75,75),DV(75,75),FV(75),FVP(75), 

1 FX(75),FXM(75),FY(75).FYM{75),PT(75),QT ( 75) 

LOGICAL LSOLVE,LPRINT,LBLK,LSTOP 

COMMON/INDX/NF,NFMAX,NP,NRHO,NGAM,L1,L2,L3,M1,M2,M3, 

IIST.JST,ITER,LAST,TITLE(13),RELAX(13),TIME,DT,XL,YL, 

21PREF,JPREF,LSOLVE(10),LPR1NT(13),LBLK(11),MODE,NTIMES(10), 
3RHOCON,ZERO 
COMMON/CNTL/LSTOP 
COMMON/SORC/SMAX,SSUM 
COMMON/COEF/FLOW,DIFF,ACOF,ALFA 
DIMENSION U(75,75),V(75,75),PC(75,75) 

EQUIVALENCE(F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3), 
1PC(1,1)) 

DIMENSION TH(75),THU(75),THDIF(75),THCU(75),THCVS(75) 
EQUIVALENCE(X(1) ,TH(1) ) ,(XU(1),THU(1) ),{XDIF( 1) ,THDIF(1) ) , 

1(XCV(1),THCV),(XCVS(1),THCVS(1)),(XL,THL) 

DIMENSION T( 75,75),STO(75,75),I BODY( 75,75) 

EQUIVALENCE (T(1,1),F(1,1,4)) 

LOGICAL LREAD,LWRITE 

DIMENSION D(10),H(10),W(10),Q(10) 

DIMENSION TTOP(75),TSIDES(2,75),TBOT(75),QTOP(75),QSIDES { 2,75) , 
1QB0T(75),QHTR(75),THBSID(75),THTOP(75),THBOT(75) 

CHARACTER*! AREAD 

(;**********************•*******************»**** 

ENTRY GRID 

^******<^******^^*I^**1<^******★********♦******^•t***** 

write(36,*) ' ks - 10 ,ra=10e6 ' 

WRITE(36,* ) 

C SET UP THE GRID 

NlX-11 
N2X-6 
N3X-9 
N4X-11 
N5X-17 
N6X-8 
N7X-10 
NlY-8 
N2Y-12 
N3Y-15 
N4Y-15 
N5Y-11 
N6y-12 

C HEATER DIMENSIONS 
WIDTH-.02 
HI-.0078 
DEPTH-.5*H1 

C SET UP THE DIMENSIONS OF THE ENCLOSURE 

AL1-3*H1/H1 
AL2-H1/H1 
AL3-(H1-DEPTH)/H1 
AL4-DEPTH/H1 
AL5-2*H1/H1 
AL6-2*H1/H1 
AL7-H1/H1 
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BL1-2*H1/H1 

BL2-2*H1/H1 

BL3-H1/H1 

BL4-2*H1/H1 

BL5-2*H1/H1 

BL6-H1/H1 

DXl-ALl/NlX 

DX2-AL2/N2X 

DX3-AL3/N3X 

DX4-AL4/N4X 

DX5-AL5/N5X 

DX6-AL6/N6X 

DX7-AL7/N7X 

DYl-BLl/NlY 

DY2-BL2/N2Y 

DY3-BL3/N3Y 

DY4-BL4/N4Y 

DY5-BL5/N5Y 

DY6-BL6/N6Y 

C MAKING OF THE GRID 

XU(2)-0 
first-3 
LLAST-3-1+N1X 
DO 11 I-FIRST,LEAST 
XU(I)-XU(I-l)+DXl 

11 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N2X 

DO 12 I-FIRST,LEAST 
XU(I)-XU(I-l)+DX2 

12 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N3X 

DO 13 I-FIRST,LEAST 
XU(I)-XU(I-l)+DX3 

13 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N4X 

DO 14 I-FlRST,LEAST 
XU(I)-XU(I-l)+DX4 

14 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N5X 

DO 15 I-FIRST,LEAST 
XU(I)-XU(I-l)+DX5 

15 CONTINUE 
FIRST-LLAST+l 
LLAST-LLAST+N6X 

DO 16 I-FIRST,LEAST 
XU(I)-XU(I-l)+DX6 

16 CONTINUE 
FIRST-LLAST+1 
LEAST-LLAST+N7X 

DO 17 I-FIRST,LEAST 
XU(I)-XU(I-l)+DX7 

17 CONTINUE 
Ll-LLAST 
XL-XU(LEAST) 

YV(2 ) - 0 
FIBST-3 
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LLAST-FIBST-l+NlY 
DO 19 I-FIRST,LLAST 
YV(I)-YV(I-l)+DYl 

19 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N2Y 

DO 20 I-FIRST,LEAST 
YV( I )-YV( I-l)+DY2 

20 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N3Y 

DO 21 I-FIRST,LEAST 
YV(I)-YV(I-l)+DY3 

21 CONTINUE 

FIRST-LLAST+1 
LLAST-LLAST+N4Y 
DO 22 I-FIRST,LEAST 
YV(I)-YV(I-l)+DY4 

22 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N5Y 

DO 23 I-FIRST,LEAST 
YV(I)-YV(I-l)+DY5 

23 CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N6Y 

DO 24 I-FIRST,LEAST 
YV(I)-YV(I-l)+DY6 

24 CONTINUE 
Ml-LLAST 
YL-YV(LEAST) 

RETURN 

C 

ENTRY START 

C CIRCUIT BOARD CORNERS 
NA-N1X+N2X+2 
NB-N1X+N2X+N3X+N4X+1 
NC-NlY+2 

ND-N1Y+N2Y+N3Y+N4Y+1 
C TEMPERATURE PRINTOUT LOCATION 
ITL-INT(N3Y/2)+l 
C HEATER LOCATIONS 
NQX-NB 
NQY1-NC+N2Y 

C INITIAL CONDITIONS 

DO 100 J-1,M1 
DO 100 I-l,LI 
U(I,J)-0. 

V(I,J)-0. 

T(I,J)-0. 

STO(I,J)»0. 

100 CONTINUE 

‘‘♦RELAXATION FOR BOUYANT FORCE (NOW FOUND AT THE TOP) 
ALFA-.1 
ALFAM--1 . -ALFA 
C PROPERTY ENTRY POINT 

RA=1.E6 
PR-25 
CP-4185 
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RHOF-999.0 
BETA-162.020E-6 
G-9.807 
CONFLD-.596 
CONPLX-10*CONFLD 
CONHTR-25*CONFLD 
DYNVIS-PR*CONFLD/CP 
XKK-CONFLD 

ALPHA-CONFLD/(RHOF*CP) 

Q1-RA*ALPHA*C0NFLD*DYNVIS/(RH0F*G*BETA*H1*H1*H1) 
USTAR-SQRT(G*BETA*Ql*Hl/CONFLD) 

AMU-(DYNVIS/(RH0F*H1*USTAR)) 

DUMMY-CP*RH0F*H1*USTAR 
C DUMMYl USED FOR DIMENSIONAL HEAT FLUXES 
DUMMY1-Q1*DUMMY/(H1*C0NFLD) 

AMUP-CONFLD/DUMMY 
GAM2-CONPLX/DUMMY 
GAM4-CONHTR/DUMMY 
XQ-01/{DEPTH*H1) 

EIN-Ql 

SOURCE-CONFLD/(CP*DEPTH*RHOF*USTAR) 

PRINT *, '1 IF WANT TO READ FILE, 2 IF STARTING FROM SCRATCH’ 

READ *, IFREAD 
IFREAD-2 

IF(IFREAD.EQ.1)THEN 
PRINT*, 'READING DATA' 


***read data*** 

NTIME-TOTAL NUMBER OF TIME STEPS 
ITIME-CURRENT TIMESTEP 
rewind(46 ) 

NTIME-1 

read(46) NTIME,Ml,Ll 

PRINT DATA FOR GRID 
read(46) (YV(J),J-2,Ml) 
read(46) (XU(I),I-2,LI) 
read(46) (Y(J),J-1,Ml) 
read(46) (X(I),I-l,Ll) 

READ DATA FOR DEPENDENDT VARIABLES 
ITIME-1 

read(46) ITIME 
DO 3000 I-l,Ml 

3000 read(46) (V(J,I),J-1,Ll ) 

DO 3020 I-l,Ml 

3020 read(46) (U(J,I),J-1,Ll) 

DO 3030 I-l,Ml 

3030 read(46) (T(J,I),J-1,Ll) 

DO 3040 I-l,Ml 

3040 read(46) (P(J,I),J-1,Ll) 

DO 3050 I-l,Ml 

3050 read(46) (STO(J,I),J-1,Ll) 

DO 3060 I-l,Ml 

3060 read(46) (F(J,I,3),J=1.Ll) 

DO 3070 I-l,Ml 

3070 read(46) (I BODY(J,I),J-1,Ll ) 
close(46 ) 

ELSE 

ENDIF 
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RETURN 


ENTRY DENSE 
RETURN 

ENTRY BOUND 
DO 300 1-2,L2 
T(1,M1)-0 
T(I,l)-0 

300 CONTINUE 

DO 301 1-2,M2 
T(Ll,I)-0 
T(l,I)-0 

301 CONTINUE 
RETURN 

ENTRY OUTPUT 

C IF(ITER.NE.0)GO TO 400 

C PRINT 401 

C 401 FORMAT(' ITER',6X,'SMAX*,8X,'SSUM 7X,'V(FRONT)', 

C 1 6X,'T(HEATER)') 

C 400 PRINT 403,ITER,SMAX,SSUM,V(NQX+1,NQYl+2),T(NQX,NQYl+2) 

C 403 FORMAT<16,1P5E12.3) 

IF(ITER.NE.0)00 TO 400 
write( 36,401 ) 

401 FORMAT(' ITER',6X,'SMAX',8X,'SSUM ',7X,’V(FRONT)’, 

1 6X, 'T( THERMO ' ) 

400 IF(INT(ITER/25)*25.NE.ITER)GOTO 404 

write(36,403)ITER,SMAX,SSUM,V(NQX+1,NQYl+lTL),T(NQX,NQYl+ITL) 

403 F0RMAT(I6,1P5E12.3) 

404 CONTINUE 

1F(ITER.LT.LAST) RETURN 

CALL PRINT 

SSUM-0. 

NCOUNT-0 

982 CONTINUE 
S2-(YV(3)-YV(2))*H1 

lF(NCOUNT.EQ.1)S2-(YV(Ml)-YV(M2))*Hl 
LO 983 I-2,L1-1 

S1-(XU{I+l)-XU(I))*H1 

POWER-2*XKK*Sl*(T(I,2)*Ql/CONFLD)/S2 

IF(NCOUNT.EQ.1)POWER-2*XKK*Sl*(T(I,M2)»Ql/CONFLD)/S2 
SSUM-POWER+SSUM 

983 CONTINUE 
NCOUNT-NCOUNT+1 
IF(NC0UNT.LT.1.5) go to 982 
NCOUNT-0 

984 CONTINUE 

S2-(XU(3)-XU(2))*Hl 

IF(NCOUNT.EQ.1)S2-(XU(L1)-XU(L2))*Hl 
DO 985 I-2,M1-1 

Sl-(YV(I+l)-YV(I))*H1 

POWER-2*XKK*Sl*(T(2,I)*Ql/CONFLD)/S2 

IFfNCOUNT.EO.1) POWER-2*XKK*Sl*(T(L2,I)*01/CONFLD)/S2 
SSUM-POWER+SSUM 

985 CONTINUE 
NCOUNT-NCOUNT+1 

IF(NCOUNT.lt. 1.5) GOTO 984 
write(36 , * ) 

write(36,*) 'ENERGY IN - ',EIN 
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ENERGY OUT 


SSUN 


write(36,*) 
write(36 , * ) ' 
write(36,* ) 
write(36,*) 'SURFACE TEMPERATURES' 
writej 36,*) 

DO 987 1-1,75 
TTOP(I)-0 
TBOTd )-0 
QTOP(I)-0 
QBOT(1)-0 
QHTB(I)-0 
THTOP(I)»0 
THBOT(I)«0 
THBS1D(1)-0 
DO 987 J-1,75 

TSIDES(1,75)-0 
TSIDES(2,75)-0 
QSIDES(1,75)-0 
QSIDES(2,75)-0 

987 CONTINUE 
C 

* CALCULATION OF SUBSTRATE SURFACE TEMPERATURES 
DO 988 I-NA,NB 

TTOP(I)-(T(I,ND)*.5*YCV{ND+1)/CONFLD+T(I,ND+1)*.5*YCV(ND)/CONPLX)/ 
1(.5*(YCV(ND+1)/CONFLD+yCV(ND)/CONPLX)) 

TBOT{I)-(T(I,NC)*.5*YCV(NC-1)/CONFLD+T(I,NC-1)*.5 *YCV(NC)/CONPLX)/ 
1(.5*{YCV(NC-1)/C0NFLD+YCV(NC)/C0NPLX)) 

988 CONTINUE 

DO 989 I-ND,NC,-1 

TSIDESd,I)-(T(NA,I)*,5*XCV(NA-1)/CONFLD+T(NA-1,I)*.5*XCV(NA) 
1/CONPLX)/( .5 *(XCV(NA-1)/CONFLD+XCV(NA)/CONPLX)) 

TSIDES(2,I)-(T(NB,I)*.5*XCV(NB+l)/CONFLD+T(NB+l,I)*.5»XCV(NB) 
1/CONPLX)/(.5*(XCV(NB+1)/CONFLD+XCV{NB)/CONPLX)) 

IF((I.LT.NQYl).OR.(I.GT.NQY1+N3Y-1))G0T0989 

TSIDES(2,1)-(T(NB,I)*.5*XCV(NB+1)/C0NPLD+T(NB+1,I)*.5»XCV(NB) 
l/CONHTR)/(.5*(XCV(NB+1)/CONFLD+XCV{NB)/CONHTR)) 

989 CONTINUE 

write(36,*) 'TOP' 
write(36,4000) (TTOP(I),I-NA,NB) 
write(36,*) 'BACK ',' FRONT' 

DO 990 I-ND,NC,-1 

990 write(36,4001) TSIDES(1,I),TSIDES(2,I) 
write(36,*) 'BOTTOM' 
write(36,4000) (TBOTd ), I-NA,NB) 

4000 F0RMAT(1X,6E11.4 / 1X,6E11.4) 

4001 FORMAT(1X,E11.4,20X,E11.4) 

C 

C SUBSTRATE SURFACE HEAT FLUXES 
wr ite{36,* ) 

write(36,*) 'SURFACE HEAT FLUXES-NONDIMENSIONAL' 
write(36 , * ) 

DO 994 I-NA,NB 

QTOPd)-(T(I,ND)-T(I,ND+1))/(.5*YCV(ND)/GAM2+.5*YCV(ND+1) 

1/AMUP) 

QBOT(I)-(Td,NC)-T(I,NC-1))/(.5*YCV(NC)/GAM2+.5 *YCV(NC-1) 

1/AMUP) 

994 CONTINUE 

write(36,4000) (QTOPd),1-NA,NB) 

DO 995 I-ND,NC,-1 

QSIDES(1,I)-(T(NA,I)-T(NA-l,!))/(.5*XCV(NA)/GAM2+.5*XCV(NA-1) 
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1/AMUP) 

QSIDES(2,I)-(T(NB,I)-T(NB+l,1))/(.5*XCV(NB)/GAM2+.5*XCV(NB4l) 
1/AMUP) 

IF( ( :.LT.NQYl).OR.(I.GT.NQYl♦NBY-1))G0T0995 

QSIDES(2,1)-(T(NB,l)-T(NB+l,I))/(.5*XCV(NB)/GAM4+.5*XCV(NB+1) 
1/AMUP) 

995 CONTINUE 

DO 996 I-ND,NC,-1 

996 write(36,4001) QS1DES(1,I),QSIDES(2,I) 
write(36,4000) (QBOTtI).I-NA,NB) 

write(36,*) 

write(36,*) 'DIMENSIONAL SUBSTRATE HEAT FLUXES' 

write(36,*) 

write(36,4000) (QTOP(I)*DUMMYl,I-NA,NB) 

DO 1050 I-ND,NC,-1 

1050 write(36,4001) QSIDES(1,I)“DUMMYl,QSIDES(2,I)‘DUMMYl 
write(36,4000) (QBOT(I)*DUMMYl,I-NA,NB) 


write(36,*) 

write(36,*) 'LOCATION OF SUBSTRATE TEMPERATURES AND FLUXES' 
write(36 , * ) 

write(36,*) 'TOP ROW ','1 Y',' & X''S' 

write(36,'(1X,E11.4 / 6E11.4 / 6B11.4)') YVIND+l),(X(I),I-NA,NB) 
wr i te(36,*) 

write(36,*) 'BACK AND FRONT X''S' 

write(36,'(1X,E11.4,2X,Ell.4)') XU(NA),XU{NB4l) 

write(36,*) 'Y POSITIONS' 

DO 1020 I-ND,NC,-1 

1020 write(36,'(5X,E11.4)') Y(I) 
write(36,*) 'BOTTOM ROW Y AND X"S' 

write(36,'(1X,E11.4 / 6E11.4 / 6E11.4)') yv(NC),(X(I),I-NA,NB) 
write(36,* ) 

write(36,*) 'HEATER SURFACE TEMPERATURES' 

DO 1021 1-NA+N3X,NB 

THTOP(I)-(T(I,NQY1+N3Y)*.5*YCV(NOYl)/CONHTR+T(I,N0Yl+N3Y-l)*.5* 
1YCV(NQY1+N3Y)/CONPLX)/(.5*(yCV(N(}Yl)/CONHTB+YCV(NQYl+N3Y)/CONPLX)) 
THBOT(I)-(T(I,NQYl)*.5*YCV(NQYl-1)/CONPLX+T(I,NQY1-1)*.5* 

1YCV(NQY1)/CONHTR)/(.5*(YCV(NQYl-l)/CONPLX+YCV(NQYl)/CONHTR)) 

1021 CONTINUE 

DO 1022 I-ND,NC,-1 

THBSID(I)-(T(NA+N3X,I)*.5*XCV(NB-N4X)/CONPLX+T(NB-N4X,I)*.5 
1*XCV(NA+N3X)/CONHTR)/(.5*(XCV(NB-N4X)/CONPLX+XCV(NA+N3X)/CONHTR)) 

1022 CONTINUE 
write(36,*) 'HEATER TOP' 

write(36,4000) (THTOP(I),I-NA+N3X,NB) 
write(36,*) 'BACK FRONT' 

DO 1025 I-NOY1+N3Y-1,N0Y1,-l 
1025 write(36,4001) THBSID(I),TSIDES(2,I) 
write(36,*) 'BOTTOM' 

write ( 36,4 000) (THBOTd ), I-NA+N3X, NB ) 

HEATER SUFACE HEAT FLUXES 
write(36, * ) 

write(36,*) 'HEATER SURFACE HEAT FLUXES - NONDIMENS IONAL' 
write(36,* ) 

IHTRC-1 

DO 997 1-NOX-N4X+1,NOX 

QHTR(IHTRC)-(T(I,N0Yl4N3Y-l)-T(1,NQYl+N3Y))/(.5*YCV(NQY1+N3Y-1) 
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1 /GAM4+.5*YCV(NQY1+N3V)/GAM2) 

IHTRC-lHTRC+1 

997 CONTINUE 

DO 998 I-NQX-N4X+1,NQX 

QHTR(IHTRC)«(T(I.NQYl)-T(I,NQY1-1 ) )/( .5*YCV(NQYl) 

1 /GAM4+.5*YCV(NQY1-1)/GAM2) 

IHTRC-IHTRC+1 

998 CONTINUE 

DO 999 I-NQY1+N3Y-1,NQY1,-l 

QHTR(IHTRC)-(T(N0X'N4X+1,I)-T(NQX-N4X,I))/(.5*XCV(NQX-N4X+1) 

1 /GAM4+.5*XCV(NQX-N4X)/GAM2) 

IHTRC«IHTRC+1 

999 CONTINUE 

write(36.4000) (QHTR(I),I-l.N4X) 

HFTOT-0 

DO 1001 1«1,N3Y 

write(36,4002 ) OHTR(2*N4X+I),QSIDES(2,NQY1+N3Y-I) 
HFTOT»HFTOT+QSlDES(2,NQY1+N3Y-I)*YCV(NQY1+N3Y-I) 

1001 CONTINUE 

write(36,4000) (QHTR(N4X+I),I-1,N4X) 

4002 FORMAT(1X.E11.4,10X,E11.4) 
write{36,* ) 

wtite(36,*) 'HEAT GOING DIRECTLY TO FLUID - ',HFTOT*DUMMYl•Hi 
write(36,* ) 
write(36,* ) 

write(36,*) 'DIMENSIONAL HEATER HEAT FLUXES' 
write(36 , * ) 

write( 36,4000 ) (QHTR(I)*DUMMYl,I-l,N4X) 

DO 1051 I-1,N3Y 

1051 write(36,4002) QHTR(2*N4X+I)*DUMMY1,QSIDES(2,NQY1+N3Y-I)*DUMMYl 
write(36,4000) (QHTR(N4X+I)‘DUMMYl,I-l,N4X) 
write(36 , * ) 

write(36,*) 'LOCATION OF HEATER TEMPERATURES AND FLUXES' 
write(36 , * ) 

write(36,*) 'TOP ROW ','1 Y',' & X"S' 

write(36,'(IX,Ell.4/6E11.4/6E11.4)')YV(NQY1+N3Y), 

1(X(I),I-NA+N3X,NB) 
write(36 , * ) 

write(36,*) 'BACK AND FRONT X"S' 

write(36,'(IX,Ell.4,2X,E11.4)') XU(NA+N3X).XU(NB+1) 
write(36,*) 'Y POSITIONS' 

DO 1030 1-NQY1+N3Y-1,NQY1,-l 
1030 write(36,'(5X,E11.4)') Y(I) 

write(36,») 'BOTTOM ROW Y AND X''S' 

write(36,'(iX,Ell.4/6E11.4/6E11.4)') yV(NQYl),(X(I),I.NA+N3X,NB) 
write(3i,* ) 

HPTOP-0 

HPBOT-0 

DO 1002 I-1,N4X 

HPTOP-HPTOP+QHTR(I)*XCV(NQX-N4X+I) 

HPBOT-HPBOT+QHTR(N4X+I)*XCV(N0X-N4X+I) 

1002 CONTINUE 
HPBACK-0 

DO 1003 I-2*N4X+1,2*N4X-*-N3Y 

HPBACK-HPBACK+OHTR(I)*YCV(NQYl) 

1003 CONTINUE 

HTRTOT-HFTOT+HPTOP+HPBOT+HPBACK 

PHTRF-100*HFTOT/HTRTOT 

PHTRT-100*HPTOP/HTRTOT 

PHTRL-100*HPBACK/HTRTOT 
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PHTRB-100*HPBOT/HTRTOT 

PHTOT-PHTRF+PHTRT+PHTRL+PHTRB 


MAXIMUM TEMPERATURE 


vir i te ( 36 , * ) 

write(36,*) 'HEATER TEMPERATURES' 
write(36 , * ) 

TMAX-0 

IXX-0 

JYY-0 

DO 1005 J-N3Y,1,-1 

wtite(36,4000) (T(NQX-N4X+I,NQY1-1+J),I-l,N4X) 
DO 1005 K-1,N4X 

IF{TMAX.LT.T(NQX-N4X+K,NQY1-1+J)) THEN 
TMAX-T{NOX-N4X+K,NOYl-l+J) 

IXX-K 
JYY-J 
ELSE 
ENDIF 
1005 CONTINUE 


1010 

1012 

1014 

1016 


write(36,* ) 

write(36, ' (1X,A,3X,E11.4)' ) 'MAX TEMPERATURE IS',TMAX 
write(36,• ) 

write(36,*) 'POSITION FROM LOWER LEFT IS' 
write(36,*) 'RIGHT',IXX 
write(36,*) 'UP',JYY 
write(36,*) 


write(36,* ) 
write(36,*) 
write(36,* ) 
write(36,* ) 
write(36 , * ) 
write(36, * ) 
write(36,* ) 
write(36 , * ) 
write(36 , * ) 
write(36,*) 


'PERCENT HEATER ENERGY LOST TO FLUID - ',PHTRF 

'PERCENT HEATER ENERGY LOST TO TOP - ',PHTRT 

'PERCENT HEATER ENERGY LOST TO BACK - ',PHTRL 

'PERCENT HEATER ENERGY LOST TO BOTTOM « ',PHTRB 


' TOTAL PERCENTAGE - ',PHTOT 


'ENERGY BALANCE' 


OLEFT-O 

DO 1010 1-2,M2 

QLEFT-QLEFT+T(2,I)*Q1*YCV(I)/XDIFl2) 

QUPPER-0 

DO 1012 1-2,L2 

(JUPPER-OUPPER+T( I ,M2)*Q1*XCV( I )/YDIF(Ml ) 

QRIGHT-0 

DO 1014 1-2,M2 

QRIGHT-QRIGHT+T(L2,I)*Q1*YCV(I)/XDIF(LI) 

QLOWER-0 

DO 1016 1-2,L2 

QLOWER-QLOWER-t-T( I ,2 ) *Q1*XCV( I )/YDlF( 2 ) 

QOUT-QLEFT+QUPPER+QRIGHT+QLOWER 

PLEFT-QLEFT*100/QOUT 

PUPPER-QUPPER*100/QOUT 

PRlGHT-QRIGHT*100/0OUT 

PLOWER-QLOWER*1OO/OOUT 

PTOT-PLEFT+PUPPER+PRIGHT+PLOWER 


write(36 , * ) 
write(36 , * ) 
write(36 , * ) 
write(36 , * ) 
write(36 , * ) 


'ENERGY IN EQUALS',EIN 
'ENERGY OUT EQUALS',QOUT 
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write(36, *) 

''“^*<36,*) 'PERCENTAGE OF HEAT LOSS THRU ENCLOSURE WALLS' 
wr 1 te(36, * ) 

write(36,*) 'LEFT WALL'.PLEFT 
write(36,*) 'TOP WALL',PUPPER 
write(36,*) 'RIGHT WALL',PRIGHT 
write(36,*) 'BOTTOM WALL',PLOWER 

write(36,*) '-- 

write(36,*) 'TOTAL EQUALS',PTOT 
write(36 , * ) 

write(36,*) 'PARAMETERS' 
write(36,* ) 

write(36,*) 'H/Hl - ',BL1+BL2+BL3+BL4+BL5+BL6 
write(36,*) 'D/Hl - ',bl 2 
write(36,*) 'L/Hl - ',BL4 
write(36,*) 'ST/Hl - ',BL5+BL6 
write(36,*) 'SB/HI - ',BLl 
write(36,*) 'DT/Hl - ',AL3+AL4 
wr i te ( 36,*) 

write(36, * ) 'RA - ',ra 
write(36,*) 'PR - ',pr 
write(36,*) 'RS - ',CONPLX/CONFLD 
write(36,*) 'RC = ',CONHTR/CONFLD 
write(36,*) 'W/Hl - ',DEPTH/HI 
write(36,*) 

***PREPARE PLOTTING DATA*** 

NTIME-TOTAL NUMBER OF TIME STEPS 
ITIME-CURRENT TIMESTEP 
NTIME-1 

write(46) NTIME,M1,L1 

PRINT DATA FOR GRID 
write(46) (YV(J),j-2,Ml) 
write(46) (XU(I),1-2,LI) 
write(46) (Y(J),J-1,Ml) 
write(46) (X(1),I-1 ,l1) 

PRINT DATA FOR DEPENDENDT VARIABLES 
ITIME-1 

write(46) ITIME 
DO 3001 I-1,M1 
3001 write(46) (V(J,I),j-1,Ll) 

DO 3021 I-1,M1 

3021 write(46) (U(J,I),j-1,Ll) 

DO 3031 I-1,M1 

3031 write(46) (T(0,I),j-1,Ll) 

DO 3041 I-1,M1 

3041 write(46) (P(0,I),j-1,Ll) 

DO 3051 I-1,M1 

3051 write(46) (STO(J,l),j-l,Ll) 

DO 3061 l-l,Ml 

3061 write(46) (F(J,I,3),J-1,Ll) 

DO 3071 I-1,M1 

3071 write(46) (IBODY(J,I),J-1,Ll) 
close(46) 

C 

RETURN 

C 

ENTRY GAMSOR 
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1F(NF.EQ.3) return 
DO 500 

DO 500 I-1,L1 
GAM(I,J)-AMU 
IBODY(1.J)-l 
IF(NF.NE.4) GO TO 500 
GAM(I,J)-AMUR 

500 CONTINUE 

DO 501 X-NA,NB 
DO 501 J-NC,ND 
GAM(I,J)-1.0E5 
IBODy(I,J)-0 
IF(NF.NE.4) GO TO 501 
GAM(I,J)-GAM2 

501 CONTINUE 
IF(NF.NE.2) GO TO 503 
DO 502 J-2,M2 

DO 502 1-2,L2 

IF(J.EQ.2)GO TO 502 

TM-(FY(J)*T(I,J)+FyM(J)*T(I,J-1)) 

CON(I,J)-ALFA*TM+ALFAM*STO(I,J) 

STO(I,J)-TM 

502 CONTINUE 

503 CONTINUE 
IF(NF.EQ.4) THEN 
DO 505 J-1,N4X 
DO 505 I-0,N3y-l 

GAM(NQX-N4X+J,NQY1+I)-GAM4 
C0N(NQX-N4X+J,NQY1+I)-S0URCE 
505 CONTINUE 
ELSE 
ENDIF 
RETURN 
END 
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APPENDIX B 

CONTOUR PROGRAM 


xxxxxxxxx PLOTTING ROUTINE FOR A SERIES OF CONTOUR PLOTS ****** 

XXKXXXKXXKXXKXXXKXXKKXXKXKXKKXMKKKXXXXKKXKKKKKXKXKXKXXKKKXXKXXXKX 

K .... CURRENT NUMBER OF TIME STEP BEING READ 

NT - TOTAL NUMBER OF TIME-STEPS ON DATA FILE 

NPTS ... TOTAL NUMBER OF POINTS (INPUT) 

NTYP ... CURVE-TYPE 

1 ... PRESSURE 

2 ... STREAMFUNCTION 

3 ... U-VELOCITY 

4 ... V-VELOCITY 

5 ... PRESSURE 


COMMON /INDATA/ X{102),Y(102),Z(102,102) 

COMMON /GRIDII/ XX(200),YY(200) 

COMMON /PLOTDT/ XPAGE,YPAGE.XPLOT,YPlOT 

COMMON /BODCTR/ XCT<102),YCT(102),XBL(102),YBL(102),NCT 

COMMON HRK (20000) 

DIMENSION XI(10404).Y1(10404),Z1(10404),ZZ(200,200) 
DIMENSION U(102,102),V(102,102),T(102.102),P{102,102), 

2 1BODY(102,102),ST1(102,102),ST2(102,102) 

DIMENSION XU(101),YV(101) 

REALX8 XARRAY(102) 

CHARACTERXl IDECl 

LOGICAL ZERO,CONT.ZERI,CONTI 


DEFINE SOME FUNCTIONS 

2ERO(II,JJ) = ABS(U(II,JJ)).LT.0.1E-]0 .AND. 

2 ABS(V(II,JJ)).LT.0.1E-10 

C0NT(I1,JJ) = (ZERO(Il+l,JJ+l).AND.ZERO(n+l,JJ-l).AND. 

2 ZEROdl-l, JJ + 1) .AND.ZEROCII-l , JJ-D) 

3 .NEQV. (ZERO(II+1,JJ+1).OR.ZERO(II+1,JJ-1).OR. 

4 ZER0(1I-1,JJ+1) .OR.ZEROdl-l, JJ-D) 
ZERIdl.JJ) = IBODYdI,JJ).EQ.O 

CONTIdDJJ) * (ZERIdl + l,JJ + 1).AND.ZERIdI + l,JJ-l) .AND. 

2 ZERIdI-1, JJ+1) .AIiD.ZERIdI-1, JJ-D) 

3 .HEQV. (ZER1(II + 1,JJ + D .OR.ZERIdI + 1, JJ-D .OR. 

4 ZERIdI-1, JJ + D.OR.ZERIdI-1, JJ-D) 

XXX DEFINE PAGESIZE AND PLOTSIZE *xx 

XPAGE =9.5 
YPAGE =11.0 
XPLOT = 6.5 
YPLOT = 6.5 

XXX read data FOR GRID xxx 
870 REWIND (3) 

READ (3) NT,LI,Ml 
READ (3) {XARRAY(D,I=2,LD 
DO 8002 1=2,LI 
002 X'Jd) = SNGL(XARRAY(D) 

READ (3) (XARRAY(J),J=2,MD 
DO 8004 J=2,M1 

8004 YV(J) = SNGL(XARRAY(J)) 

READ (3) (XARRAYd),I = l,LD 
DO 8006 1=1,LI 
8006 Xd) = SNGL(XARRAY(D) 

READ (3) (XARRAY(J),J=I,MD 
DO 8008 J=1,M1 
8008 Y(J) = SNGL(XARRAY(J)) 

PREPARE THE GRID-DATA FOR THE CONTOUR PLOTTING 
NPTS = LlxMl 
DO 1300 J=1,M1 
DO 1300 1=1,LI 
Jl=Llx(J-l)+I 
1300 XKJD = Xd) 

DO 1310 J=1,M1 
DO 1310 1=1,LI 
Jl = Llx( J-D + I 
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1310 Yl(Jl) = Y(J) 

C ««* NOMINATE DEVICE «»» 

C 

WRITE (>(,930) 

930 FORMAT (IHl,3A(1H»)/lX,’PLEASE, SELECT PLOTTING DEVICE'/IX,34(IK*: 


2 //5X,'TEKTRONIX 618 SCREEN ... (DV 

3 5X,'PRINTABLE SHERPA FILE ... (2)'/ 

4 SX.'IBM COLOR MONITOR . (3)'/ 

6 5X.'COMPRESSED METAFILE . (4)'/ 

5 SX.'EXIT PLOTTING . (99)') 

READ (K,*) IPLOT 

IF ((IPLOT.LT.n.OR (IPLOT.GT.4)) GOTO 180 
IF (IPLOT.EQ.l) CALL TEK618 


IF (IPLOT.£0.2) CALL SHERPA (’CONTOUR ’,’T’,2) 

IF (IPLOT.EQ.3) CALL IBM79 
IF (IPLOT.EO.4) CALL COMPRS 
C 

C »** GET THE NUMBER OF TIME-STEPS TO BE PLOTTED AND THE INCREMENT *** 
WRITE (>(,910) NT 
910 FORMAT (1 HI,39(1H*)/, 

2 IX,'YOUR INPUT FILE INCLUDES ',13,' TIME-STEPS'/IX.39(IHK)// 

3 IX,'CHOOSE A METHOD FOR CONSTRUCTING A REGULAR MATRIX;'/ 


4 5X,'UNIFORM GRID/NO REFINEMENT (102X102). (0)'/ 

5 5X,’DISSPLA"S GETMAT ALGORITHM ( 102x102). ( 1 )'/ 

6 5X,'CUBIC SPLINE INTERPOLATION (102x102). (2)'/ 

7 5X,'LINEAR INTERPOLATION (200X200). (3)') 


READ (x,«) IMAT 
IF (IMAT.NE.O) GOTO 202 
LL = LI 
MM = Ml 
GOTO 204 

202 WRITE (x, 9 ll ) 

911 FORMAT (IX,'ENTER THE DIMENSIONS OF THE MATRIX TO BE CONSTRUCTED;' 

2 ,' (IXDIM.IYDIM)') 

READ (X,X) LL,MM 

C GENERATE AN EQUALLY SPACED GRID 
204 DELX = X(L1)-X(1) 

DELY = Y(M1)-Y(1) 

DO 206 1 = 1 ,LL 

206 XX(I) = X(l) + DELXXFL0AT(I-1)/FL0AT(LL-1) 

DO 208 J=1,MM 

208 YY(J) = Y(l) + DELYxFL0AT(J-1)/FL0AT(MM-1) 

c_ 

C *** LOOP OVER A SERIES OF PLOTS xxx BEGIN OF "20O-LOOP" 

C 

DO 200 N=1,NT 
READ (3) K 
WRITE (X.912) K 

912 FORMAT (//IX,'DO YOU WANT ANY PLOTS FOR THE ',13, 

2 '-TH TIME-STEP’ (Y/N)') 

READ (x,914) IDECl 
914 FORMAT (Al) 

IF (IDECl.EQ.'Y') GOTO 212 
C ADVANCE BY 1 TIME-STEP 
DO 210 1=1,6 
210 READ (3) 

GOTO 200 

C READ DEPENDENT VARIABLES FROM FILE 3 

212 DO 213 1 = 1 ,LI 

READ (3) (XARRAY(J),J=1,M1) 

DO 213 J=1,M1 

213 U(I,J) = SNGL(XARRAY(J)) 

DO 214 1 = 1,LI 

READ (3) (XARRAY(J),J=1,M1) 

DO 214 J=1,M1 

214 V(I,J) = 5NGL(XARRAY(J)) 

DO 215 1 = 1,LI 

READ (3) (XARRAY(J),J=1,H1) 

DO 215 J=1,M1 

215 T(I,J) = SNGl(XARRAY(J)) 

DO 216 1 = 1,LI 

READ (3) (XARRAY(J),J=1,M1) 
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DO 216 J=1,M1 

216 Pd.Jl = SNGKXARRAYCJ)) 

DO 217 1 = 1,LI 

READ (3) (XARRAYCJ),J=1,M1) 
DO 217 J=1,M1 

217 STICI.J) = SNGKXARRAYCJ)) 
DO 218 1=1,LI 


READ C3) (XARRAY(J),J=1,M1) 

DO 218 J=1,M1 

218 ST2(I,J) = SNGKXARRAYCJ)) 

DO 219 1=1,LI 

219 READ C3) C1BODYC1,J),J=1,Ml) 

IF CN.NE.l) GOTO 890 
WRITE CX,939) 

939 FORMAT C/IX,'SHOULD THE CONTOUR OF THE BODY BE IDENTIFIED VIA:'/ 

2 6X,'THE VELOCITIES . Cl)'/ 

3 6X, 'THE INPUT DATA . C2) ' ) 

READ CK,*) ICONT 

IF CIC0NT.E(}.2) GOTO 1500 


xxx FIND THE CONTOUR OF THE BODY VIA THE VELOCITIES 

LOCATE A CONVEX CORNER OF THE BODY 
DO 1320 J=A,Ml-3 
DO 1320 I=4,Ll-3 
IF C2ER0CI,J)) GOTO 1322 
1320 CONTINUE 
1322 ICT = 1 

XCTCICT) = XCI) 

YCTCICT) = YCJ) 

IP = 1 

JP = 0 

DO 1390 ICT=2,102 
DO 1329 IDUM=1,MAXCL1,M1) 

IF CZEROCI+IP,J+JP).AND.COHTCI+IP,J+JP)) GOTO 1330 
IF CIP.EQ.O) THEN 
JP = 0 

DO 1326 IP=-1,1,2 

1326 IF CZER0(I+IP,J+JP).AND.C0NTCI+1P,J+JP)) GOTO 1332 

ELSE 

IP = 0 

DO 1328 JP=-1,1,2 

1328 IF CZEROCI+IP,J+JP).AND.CONTCI+IP,J+JP)) GOTO 1532 

ENDIF 

ANOTHER CONTOUR POINT HAS BEEN FOUND 
1330 I = I+IP 
J = J+JP 

1329 CONTINUE 

THE PREVIOUSLY FOUND POINT IS A CORNER 
1332 XCTCICT) = XCI) 

YCTCICT) = YCJ) 

IF CABSCXCTCICT)-XCTC1)).LT.O.IE-6.AND. 

2 ABSCYCTCICT)-YCTCl)).LT.0.1E-6) GOTO 1350 
I = I+IP 
J = J+JP 
1390 CONTINUE 
1350 NCT = ICT 
GOTO 1900 


«*» FIND THE CONTOUR OF THE BODY VIA INPUT *** 

LOCATE A CONVEX CORNER OF THE BODY 
1500 DO 1520 J=9,Ml-3 
DO 1520 I=9,Ll-3 
IF C2ER1CI,J)) GOTO 1522 
1520 CONTINUE 
1522 ICT = 1 

XCTCICT) = XCI) 

YCTCICT) = YCJ) 

IP = 1 
JP = 0 

DO 1590 ICT=2,102 
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GOTO 1550 


DO 152'i IDUM=:l.MAX(Ll,m) 

IF (ZERKI + IP, J + JP) .AND. CONTI (I + IP,J + JP ) ) 

IF (IP.EQ.O) THEN 
JP = 0 

DO 1526 IP=-1,1.2 
1526 IF C2ERI(I+IP,J+JP).AND.CONTKl+IP,J+JP)) GOTO 1532 
ELSE 

IP = 0 

DO 1528 JP=-1.1,2 

1528 IF (ZERKl+lP,J+JP).AND.C0NT1(1+1P,J+JP)) GOTO 1532 
ENDIF 

C ANOTHER CONTOUR POINT HAS BEEN FOUND 
1530 I = 1+IP 
J = J+JP 
1524 CONTINUE 

C THE PREVIOUSLY FOUND POINT IS A CORNER 
1532 XCT(ICT) = X(I) 

YCT(ICT) = Y<J) 

IF (ABS(XCT(ICT)-XCT(l)).LT.0.1E-6.AND. 

2 ABS(YCT(ICT)-YCT(l)).LT.0.1E-6) GOTO 1550 
1 = 1+IP 
J = J+JP 
1540 CONTINUE 
1550 NCT = ICT 

COfIPUTE THE COORDINATES FOR BLANKING (IN INCHES FROM PHYSOR) 
1400 DO 1410 ICT=1,NCT 

XBL(ICT) = (XCT(ICT)-X(1))/(X(L1)-X(1)))(XPL0T 
1410 YBL(ICT) = (YCT(1CT)-Ytl))/(Y(M1)-Y(1))*YPL0T 

xkk find out WHICH CURVES TO PLOT kkx 
890 WRITE (»,900) 

900 FORMAT (IHl.'WHAT WOULD YOU LIKE TO PLOT?'/ 


2 5X,'TEMPERATURE. (1)'/ 

3 5X,'STREAMLINES. (2)'/ 

4 5X,'U-VELOCITY . (5)'/ 

5 5X,'V-VELOCITY . (4)'/ 

6 5X,'PRESSURE . (5)'/ 

7 5X,'NEXT TIMESTEP . (0)'/ 

8 5X,'RE-INITIALISE . (-99)'/ 

9 5X,'EXIT PLOTTING . (99)' ) 

READ (K,«) NTYP 


IF (NTYP,EQ.99) GOTO 180 
IF (NTYP.EQ.-99) GOTO 870 
IF (NTYP.EQ.O) GOTO 200 

ASSIGN THE VALUES OF THE Z-FUNCTION WITH THE CHOSEN DATA 
GOTO (1100,1102,1104,1106.1100 NTYP 

TEMPERATURE 

1100 DO 1101 J=1,M1 
DO 1101 1=1,LI 

1101 Z(I,J) = T(I,J) 

GOTO 1120 

C STREAMFUNCTION 

1102 DO 1103 J=1,M1 
DO 1103 1=1,LI 

1103 Z(1,J) = ST2(I,J) + .0005222 
GOTO 1120 

C U-VELOCITY 

1104 DO 1105 J=1,M1 
DO 1105,1=1,LI 

1105 Z(I,J) = U(I,J) 

GOTO 1120 

C V-VELOCITY 

1106 DO 1107 J=1,M1 
DO 1107,1=1,LI 

1107 Z(I,J) * V(I,J) 

GOTO 1120 

C PRESSURE 

1108 DO 1109 J=1,M1 
DO 1109 1=1,LI 
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1109 Z(I.J) = Pd.J) 

COMPUTE AND DRAW CONTOUR-LINES 
1120 CALL CONTDD (NTYP,IMAT,XI.Y1,Z1,ZZ,L1,Ml,L1xMl,LL,MM) 
GOTO 890 

200 CONTINUE 

END OF "200-LOOP" 


EXIT PLOTTING ROUTINE xxx 
180 CALL DONEPL 
STOP 
END 




SUBROUTINE CONTDD (NTYP,IMAT,XI,Y1,21,ZZ,L1,Ml,NPTS,LL,MM) 

THIS SUBROUTINE COMPUTES A REGULAR MATRIX (ON AN EQUALLY SPACED 
GRID) EITHER BY DISSPLA'S GETMAT-ALGORITHM OR BY A CUBIC SPLINE 
INTERPOLATION. SUBSEQUENTLY CONTOUR-LINES ARE DRAWN. 

KlOCKXKXXXXXXXKXKXIOOiXXXXXKXKXXXXKKXXXXXXXXKXXKKXKXXKXifXXiOOOOfXXXiOOOOO 

COMMON /BODCTR/ XCT(102).YCT{102),XBL(102),YBL(102),NCT 
COMMON /INDATA/ X<102).Y(102),2<102,102) 

COMMON /GRIDII/ XX(200),YY(200) 

COMMON /PLOTDT/ XPAGE,YPAGE,XPLOT,YPLOT 
COMMON WRK (20000) 

DIMENSION XI(NPTS),Y1(NPTS),Z1(NPTS),ZZ(LL,MM) 

DIMENSION ZI(200,200),C(<i,102) 

CHARACTERXl LCONT,LGRID 


XXX INPUT THE INCREMENT BETWEEN CONTOUR LEVELS xxx 
WRITE (X,900) 

900 FORMAT (IX,'SHOULD DISSPLA PICK THE INCREMENT BETWEEN CONTOUR ', 
2 'LEVELS? (Y/N)') 

. READ (»,910) LCONT 

910 FORMAT (Al) 

IF (LCONT.EQ,'Y') GOTO 925 
WRITE(x,920) 

920 FORMAT (IX,'INPUT THE INCREMENT BETWEEN CONTOUR LEVELS:') 

READ (x,x) ZINCR 

925 WRITE (x,930) 

930 FORMAT (IX,'DO YOU WANT A GRID TO BE DRAWN? (Y/N)') 

READ (x,910) LGRID 

xxx SET UP ALPHABET xxx 
CALL RESET CALL') 

CALL SWISSM 
CALL HWSHD 
CALL CHRPAT(16) 

CALL BASAIF ('L/CSTD') 

CALL MIXALF ('STAND') 

CALL SETCLR ('BLUE') 

CALL SHDCHR(90.,1,.002,1) 

xxx DEFINE PLOT AREA xxx 

CALL PAGE (XPAGE,YPAGE) 

IF (IPLOT.EQ.l) CALL HWSCAL ('SCREEN') 

CALL PHYSOR (1.25,1.0) 

CALL AREA2D (YPLOT,XPLOT) 

xxx TEMPORARY FIX FOR CONTOURING PROBLEMS 
1 = 10 
IF 1 = 1 

CALL THKFRM (.01) 

xxx LABEL AXES xxx 

CALL HEIGHT (0.18) 

CALL SETCLR ('GREEN') 


110 






non nr> oooo oon ooo o on 


IF (NTYF.EO.I) CALL HEADIN (•(T)EMPERA1URE (C.)0NT0UR5$’, 
1 100,1.5,1) 

IF (NTYP.E0.2) CALL HEADIN (•(SITKEAML1 NEDS', 

1 100,1.5,1) 

IF (HTYP.EO.S) CALL HEADIN (•(U-V)ELOCITY (C)0NT0UR5«', 

1 100,1.5,1) 

IF (NTYP.EO.A) CALL HEADIN ( • (V-V)ELOCITY (OONTOURSi ' , 

1 100,1.5,1) 

IF (NTYP.E0.5) CALL HEADIN ('(PIRESSURE (C)0NT0UP.5$ • 

1 ,100,1.5,1) 

CALL RESET ('HEIGHT') 

CALL SETCLR ('BLUE') 


nn* SET UP AXES »»* 

CALL XNAME ('Y $' ,100) 

CALL YAXANG (0.0) 

CALL YNAME ('X $' ,100) 

XINC = 0.1«(X(L1)-X(1)) 

YIMC = 0.1»(Y(M1)-Y(1) ) 

CALL GRAF (Y(1),YINC,Y(M1),X(1),XINC,X(L1)) 

CALL GRID (0,0) 

CALL BCOMON (20000) 

DRAW A GRID (IF REQUESTED) 

IF (LGRID.EQ.'N') GOTO 90 
DO 30 1=1,LI 

50 CALL RLVEC (Y(1),X(I),Y(M1),X(I),0000) 

DO 32 J=1,M1 

32 CALL RLVEC (Y(J),X<1),Y(J),X(M1),0000) 

PREPARE CONTOURS PLOTS xx* 

90 GOTO (100,200,300) IMAT 

XXX INPUT MATRIX IS REGULAR / NO GRID REFINEMENT REQUIRED xxx 

00 10 J=1,MM 
DO 10 1=1,LL 
10 2Z(I,J) = 2(J,1) 

GOTO 900 

XXX CONSTRUCT THE MATRIX BY DISSPLA'S ALGORITHM GETMAT xxx 

TRANSFORM THE Z 2-D FUNCTION TO 1-D FUNCTION 
100 DO no J = 1,M1 
DO no 1 = 1.LI 
Jl=Llx{J-l ) + I 

no zKji) = z(i,j) 

CONSTRUCT A REGULAR MATRIX 
CALL BGNMAT(MM,LL) 

CALL GETHAT(Y1,X1,Z1,NPTS.0) 

CALL ENDMAT(ZZ,0) 

GOTO 900 

xxx CONSTRUCT THE MATRIX BY CUBIC SPLINE INTERPOLATION xxx 

200 DO 210 J=1,M1 
DO 212 1=1,LI 
212 C(1,I) = Z(I,J) 

CALL CUBSPL (X,C,L1,0,0) 

DO 219 1 = 1,LL 

219 Z1(I,J) = PPVALU (X,C,L1,9,XX(I),0) 

210 CONTINUE 

DO 220 1 = 1,LL 
DO 222 J=1,M1 
222 C(l, J) = ZKI ,J) 

CALL CUBSPL (Y,C,M1,0,0) 

DO 229 J=1,MM 

229 ZZ(J,I) = PPVALU (Y,C,M1,9,YY(J),0) 

220 CONTINUE 
GOTO 900 
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c 

C CONSTRUCT THE MATRIX BY LINEAR INTERPOLATION »»» 

C 

300 DO 310 J=1.M1 
DO 310 1 = 1 . LL 

CALL INTERV (X,L1,XXC1),LEFT,MFLAG) 

310 Z1(I,J) = ZCLEFT.J) + (Z(LEFT+l,J)-Z(LEFT,J))*(XXCI)-X(LEFT))/ 

2 <X(LEFTa)-X(LEFT)) 

C 

DO 320 1=1,LL 
DO 320 J=1,MM 

CALL INTERV (Y,M1.YYCJ),LEFT,MFLAG) 

320 ZZ(J,I) = ZKI.LEFT) + ( ZI (I, L EFT +1)-ZI (I , L EFT ) )«(YYC J )-Y( L EFT ))/ 

2 (YCLEFT + D-YCLEFT)) 

C 

C »«* DRAW THE BODY CONTOUR AND BLANK THE INSIDE *** 

AOO CALL LINEAR 

CALL SETCLR ('CYAN') 

CALL THKCRV (0.025) 

CALL CURVE (YCT,XCT,NCT,0) 

CALL RESET CTHKCRV) 

IF (NTYP.EO.l .OR. LGRID.EO.'Y') GOTO AlO 
C CALL SHADE (YCT,XCT,NCT,A5.,.1,1,0,0) 

CALL BLPOLY (YBL,XBL.NCT,1) 

C 

C K** COMPUTE THE CONTOUR-LINES K** 

AlO CALL HEIGHT (0.12) 

CALL SETCLR ('RED') 

IF (LCONT.EQ.'Y') CALL CONMAK (ZZ,MM,LL,5HSCALE) 

IF (LCONT.EO .'N') CALL CONMAK (ZZ,MM,LL,ZINCR) 

CALL COHLIN(0,5HSOLID,6HLABELS,1,5) 

C 

C »** DRAW CONTOUR LINES «*)( 

CALL C0NANG(135.) 

CALL CONMINCl.S) 

C CALL RASPLN(0.25) 

CALL C0NTUR(1,6HLABELS,AHDRAW) 

CALL RESET CTHKCRV) 

C 

C TERMINATE THE CURRENTLY DISPLAYED GRAPH 
CALL ENDGR(O) 

CALL ENOPLCO) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCt 


c ( 

SUBROUTINE CUBSPL (TAU,C.N,IBCBEG,IBCEND) 

C ( 

C THIS SUBROUTINE COMPUTES THE CUBIC SPLINE INTERPOLANT FOR ( 

C DIFFERENT BOUNDARY CONDITIONS. C 

C SOURCE: CARL DE BOOR, A PRACTICAL GUIDE TO SPLINES C 

C C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C TAU . ABSCISSAE OF THE DATA POINTS 

C C(l,I) .... ORDINATES OF THE DATA POINTS 

C N . NUMBER OF DATA POINTS 

C IBCBEG .... INDICATOR OF BOUNDARY CONDITION AT TAU(l) 

C IBCBEG=0: NOT-A-KNOT-CONDITION 

C IBCBEG=1: PRESCRIBED SLOPE 

C IBCBEG=2; PRESCRIBED SECOND DERIVATIVE 

C IBCEND .... INDICATOR OF BOUNDARY CONDITION AT TAU(N) 

C C(J,I) _ POLYNOMIAL COEFFICIENTS OF SPLINE INTERPOLANT 

C F(X) = C(l,I) + H«(C(2,I) + H)((C(3,I)-tH)(C(4,I)/3. )/2. ) 

C WHERE H = X - TAU(l) 

C _ 

C 

DIMENSION C(A,N),TAU(N) 

C _ 

C 

L = N-1 

C COMPUTE FIRST DIFFERENCES AND FIRST DIVIDED DIFFERENCES 
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C (C SERVES AS A TEMPORARY STORAGE) 

DO 10 M=2,N 

C(3,M) = TAU(M) - TAU<M-1) 

10 C(A,M) = (C(l.M) - C(1.M-1))/C(3,M) 

CONSTRUCT FIRST EQUATION FROM FIRST BOUNDARY CONDITION «»» 

IF (IBCBEG-1) 11,15.16 

11 IF (N.GT.2) GOTO 12 

NOT-A-KNOT COND AT LEFT END, AND N.EO.2 
C(A,1) =1.0 
C{3,1) = 1.0 
C(2,l) = 2.0*C(4,2) 

GOTO 25 

NOT-A-KNOT COND AT LEFT END, AND N.6T.2 

12 CC4,1) = C(3,3) 

C(3,l) = C(3,2) + C(3^3) 

C(2,l) = ((C(3,2)+2.0)«C(3.1))«C(4,2)«C(3,3)+C(3,2)*«2*C(4,3)) 

2 /C(3,l) 

GOTO 19 

SLOPE PRESCRIBED AT LEFT END 

15 C(4,l) = 1.0 
C(3,l) = 0.0 
GOTO 18 

SECOND DERIVATIVE PRESCRIBED AT LEFT END 

16 C(4,l) = 2.0 
C(3,l) = 1.0 

C(2,l) = 3.0*C(4,2) - C(3,2)/2.0»C(2.1) 

18 IF <N.EQ.2) GOTO 25 

IN CASE OF INTERIOR KNOTS SET UP EQU, EXECUTE FORWARD PASS OF GAUSS 

19 DO 20 M=2,L 
G = -C(3,M4-1)/C(4,M-1) 

C(2,M) = GKC(2,M-1) + 3.0X(C(3,M"'C ■, M+1) +C ( 3, M+1) >(C( 4, M)) 

20 C(4,M) = G*C(3,M-1) + 2.0X(C(^,M) + C(5,M*1)) 

KXX CONSTRUCT LAST EUQATION FROM SECOND BOUNDARY CONDITION 
IF (IBCEND-1) 21.30,24 

FOR PRESCRIBED SLOPE ENTER DIRECTLY THE BACK SUBSTITUTION 

21 IF ((N.E0.3) .AND.dBCEEG EQ 0)) GOTO 22 

NOT-A-KNOT COND AT RIGHT END, AND N.GE.3, 

AND (CN.GT.3) .OR. (NOT-A-KNOT COND AT LEFT END)) 

G = C(3,N-1) + C<3.N) 

C(2,N) = ((C(3,N)+2.0*G)»C{4,N)*C(3,N-1) 

2 + C<3,N)«*2*(C(l,N-l)-C(l,N-2))/C(3,N-l))/G 

G = -G/C(4,N-1) 

C(4,N) = C(3,N-1) 

GOTO 29 

NOT-A-KNOT COND AT RIGHT END, AND C(N . EQ . 3). AND.(NOT-A-KNOT AT LEFT) 
.OR. (<N.E0.2).AND.(NOT-A-KNOT COND AT LEFT END)) 

22 C(2,N) = 2.0*C(4,N) 

C(4,N) =1,0 
GOTO 28 

SECOND DERIVATIVE PRESCRIBED AT RIGHT END 

24 C(2,N) = 3.0*C(4,N) + C( 3, N)/2.0»(C( 2, N ) 

C(4,N) = 2.0 
GOTO 28 

25 IF (IBCEND-1) 26,30,24 

26 IF (IBCBEG.GT.O) GOTO 22 

NOT-A-KNOT COND AT RIGHT AND LEFT END, AND N.E0.2 
C(2,N) = C(4,N) 

GOTO 30 

28 G = -1.0/C(4,N-I) 

XXK SOLVE TRIDIAGONAL SYSTEM FOR THE SLOPES C(2,J) 
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C FORWARD PASS (=ELIMINATION OF LOWER DIAGONAL) 

29 C(A,N) = G«C(3,N-1) + C(4.N) 

C(2,N) = tG*C(2.N-l) + C(2.N))/C(4,N) 

C BACKWARD PASS (=BACK SUBSTITUTION) 

50 DO 40 J=L,I,-1 

40 C(2.J) = (C(2.J) - C<3,J)«C(2.J+l))/C(4,J) 

C _ 

C 

C GENERATE CUBIC COEFFICIENTS IN EACH INTERVAL 
C 

DO 50 1=2.N 
DTAU = C(3,I) 

DIVDFl = (C(1,I) - C(1,I-I))/DTAU 
DIVDF3 = C(2,I-1) 4 C(2,I) - 2.0XDIVDF1 
C(3.I-1) = 2.0«(DIVDF1-C(2,I-1)-D1VDF3)/DTAU 
50 CC4,I-1) = (DIVDF3/'DTAU)*(6.0/DTAU) 

C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

REAL FUNCTION PPVALU (BREAK.COEF,L.K,X,JDERIV) 

C 

C THIS FUNCTION CALCULATES THE VALUE OF THE JDERIV-TH DERIVATIVE 
C OF A PIECEWISE POLYNOMIAL FUNCTION. 

C SOURCE: CARL DE BOOR. A PRACTICAL GUIDE TO SPLINES 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c 

c BREAK . ABSCISSAE OF BREAKPOINTS 

C COEF . COEFFICIENTS OF PIECEWISE POLYNOMIAL INTERPOLANT 

C L . NUMBER OF BREAKPOINTS 

C K . ORDER OF POLYNOMIAL FUNCTIONS 

C X . ABSCISSA AT WHICH TO EVALUATE 

C JDERIV . ORDER OF DERIVATIVE (0 = FUNCTION) 

C _ 

C 


DIMENSION BREAK(L),COEF(K,L) 

C _ 

C 

PPVALU = 0.0 
FMMJDR = K - JDERIV 

C DERIVATIVES OF ORDER K OR HIGHER ARE IDENTICALLY ZERO 
IF (FMMJDR.LE.O.O) GOTO 99 
C 

IF (L.EQ.l) THEN 

C SINGLE polynomial FUNCTION 
I = 1 
ELSE 

C FIND INDEX I OF LARGEST BREAKPOINT TO THE LEFT OF X 
CALL INTERV (BREAK,L.X,I,HDUMMY) 

ENDIF 

H = X - BREAK(I) 

C 

C EVALUATE JDERIV-TH DERIVATIVE OF I-TH POLYNOMIAL PIECE AT X 
DO 10 M=K,JDERIV+L.-l 
PPVALU = (PPVALU/FMMJDR)»H ♦ COEF(M.I) 

10 FMMJDR = FMMJDR - 1.0 
C 

99 RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE INTERV (XT.LXT,X,LEFT,MFLAG) 

C 

C THIS SUBROUTINE DETERMINES THE INTERVAL INDEX WITH RESPECT 
C TO THE BREAKPOINT SEQUENCE. 

C SOURCE; CARL DE BOOR, A PRACTICAL GUIDE TO SPLINES 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c XT . ABSCISSAE OF BREAKPOINTS 
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LXT . NUMBER OF BREAKPOINTS 

X . ABSCISSA WHOSE INTERVAL INDEX IS TO BE DETERMINED 

LEFT . INDEX OF THE LARGEST BREAKPOINT TO THE LEFT OF X 

MFLAG . RETURN CODE: MFLAG =-I: X.LT.XT(l) 

MFLAG = 0: XT(l).LE.X.LT.XT(LXT) 
MFLAG = I: X.GE.XT(LXT) 


DIMENSION XT(LXT) 


RECALL THE INDEX OF THE PREVIOUS CALL TO THIS SUBROUTINE 
SAVE ILO 
DATA ILO/1/ 

IHI = ILO + I 
IF (IHI.LT.LXT) GOTO 20 
IF tX.GT.XT(LXT)) GOTO 110 
IF (LXT.LE.l) GOTO 90 
ILO = LXT -1 
IHI = LXT 

20 IF (X.GE.XT(IHI)) GOTO 40 
IF (X.GE.XTCILO)) GOTO 100 

DECREASE ILO TO CAPTURE 
ISTEP = 1 
31 IHI = ILO 

ILO = IHI - ISTEP 
IF (ILO.LE.l) GOTO 35 
IF (X.GE.XTCILO)) GOTO 50 
ISTEP = ISTEPX2 
GOTO 31 
35 ILO = 1 

IF (X.LT.XTCD) GOTO 90 
GOTO 50 

INCREASE IHI TO CAPTURE X 

40 ISTEP = 1 

41 ILO = IHI 
IHI = ILO + ISTEP 
IF (IHI.GE.LXT) GOTO 45 
IF (X.LT.XT(IHI)) GOTO 50 
ISTEP = ISTEPX2 
GOTO 41 

45 IF (X.GE.XT(LXT)) GOTO 110 
IHI = LXT 

NARROW THE INTERVAL (XT(ILO),XT(IHI)) 

50 MIDDLE = (ILO + IHI)/2 

IF (MIDDLE.EQ.ILO) GOTO 100 

IF (X.LT.XKMIDDLE)) GOTO 53 
ILO = MIDDLE 
GOTO 50 

53 IHI = MIDDLE 
GOTO 50 

SET OUTPUT AND ASSIGN RETURN CODES 

X LIES TO THE LEFT OF THE SMALLEST BREAKPOINT 
90 MFLAG = -1 
LEFT = 1 
RETURN 

C REGULAR RETURN FROM THIS SUBROUTINE 
100 MFLAG = 0 
LEFT = ILO 
RETURN 

C X LIES TO THE RIGHT OF THE LARGEST BREAKPOINT 
no MFLAG = 1 

LEFT = LXT-I 

RETURN 

END 









APPENDIX C 
TRAVERSE PROCOtAM 


10 ! !M M !!M!!!!!! !! ! ! ! ! ! ! ! ! M ! ! j ! ! ! I ! ! ! ! I ! ! ! t I ! I r ! I I M I I! ( ! ! ! ! I I M M I I 

20 ! ! PROGRAM TRAU 

30 ! ! PROGRAM USED TO MOUE TRAVERSE 

31 ! ! 

32 ! ! 

40 <!!M!I I!!I!I ! ! j !! ! !! I I ! ! H ! ! ! ! j ! I ! ! j ! M ! ! ! ! !! ! M ! j f ! ! ! ! ! !! M ! I! I ! I H 

42 Isc-9 

43 ASSIGN @Datacomnt TO Isc 

44 ASSIGN STerninal TO 1 

4G CONTROL Isc.3;1200 

47 CONTROL I sc ,4 ,• 2+4+0+1 G 

49 CONTROL Isc.5;10 

50 COM /Path/ SDataccmn 

52 OUTPUT SOatacopiMi "EZT" 

53 OUTPUT eOatacorw;"C S1000 .SM2000.11500 .IM1500 .R“ 

54 END 
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APPENDIX D 


TEMPERATURE ACQUISITION PROGRAM 

M ! ! ! M I ! I ! I ! H I I I I I M I I I I I I I I ! M M I ! M ! n I I I ! I M I I ! I I I I I I I ! I I ! ! > ! 
!! TEMPERATURE ACQUISITION PR06RAM 

!! (STEADY STATE) 

!! PROSRAM PERFORMS THE.REQUIRED FUNCTIONS NEEDED TO OBTAIN DATA 
!I FROM THE ACQUISITION SYSTEM AND CONUERT THE THERMOCOUPLE OUT- 
!! PUT UOLTASES INTO TEMPERATURE VALUES. THE TEMPERATURES ARE 
!! STORED INTO DATA FILES FOR FURTHER REDUCTION. 

! VARIABLE DEFINITIONS 
!! VOLTS : THERMOCOUPLE OUTPUT 

i! TEMP : HEATER TEMPERATURES 

M ' ‘ ! I ! > I i ! ! ! M ! ! ! M I ! I M M n ! ! ! M ! M I ' ‘ M < M ! ! I I { I ! H ! ! H M ! ! ' ' > 
pcAi_ VcIts(SO) 

REAL Tenp(£0) 

ICREATE BOAT "THRED!5: ,700,!" 4E 
! ASSIGN @PatA3 TO "THREDIE: 700 1" 

PRINT OATE$(TIMEDATE) 

PRINT " " 

PRINT “ 1.0 UATT INPUT POWER” 
p.cjTr^jT •' 45 CT.CTDC fijr, SHROUD" 

PRINT ” BOAT FILE THREDIE:,700 I" 

pOT.NjT " " 


go pOTfiJI " '• 

53 I M I ! I I I I I I I I M I ( I I I I I M M I I cyppAgg 4!-20 

9E OUTPUT 70Si "CONFMEAS DCV, 100-M9,USE 0" 

Qg pnp T.) tq 20 

9S ENTER 709!Vclt5(I) 

99 Tswp( I )• ,0005797+( 25S25,1329*Volts( I))-(607799.24E7»( Vcl 15( I )*Vcl*.: 
21 952034.3354•( Volts', I !''3 ))+(9370910995,1974*( Volt5( I )"4 ) >+ .039 

100 NEXT I 

E30 ! ! ! ! I ' I ! ! I ■ ! ! ! M ! ! I M M M ' ' SURFACE t2I-t40 

EE0 OUTPUT 709i"CONFMEAS DCV ,200-219 .USE 0” 

E50 FOR I»2! TO 40 

E70 ENTER 709iVclts(I) 

590 Te.ip( I )-.000E797+(25925,132S»VoIt5( I) )-(S07799.2467»( Volts'. I .>*'.'o:t! 

21952034,33S4*( Volts( I .'"3 )+( 9370910995.1974o( VoltsC I )'4 ) )+.039 
590 NEXT I 

500 I ! I I M I ! ! I MI I I I I I IMI I I!I SURFAC 141-145 
610 OUTPUT 709;"CONFMEAS DCV ,300-304 .USE 0" 

£20 FOR 1-41 TO 45 

955 ENTER 709iVoits(I ) 

957 Tenp'I 00067S7+(25925.1329*VoIt5(I> >-(507799.2467«(VoltsvI>»Vcits 
2!952034.3364*(Volt5(I )'3 ))+(9370910996.1974*<Volt5(I)"4 ) )+.039 
959 NEXT I 

979 I!! PRINT OUTPUT TEMPERATURES 
979 PRINT ■ " 

890 FOR 1-1 TO 45 

99! PRINT "HEATER l"il," TEMPE.RATL'RE IN DEG. C -"jTenpd) 

992 NEXT I 

993 'OUTPUT SPat.H3;Te.+p( * ) 

994 PRINT " " 

£95 PRINT " " 
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APPENDIX E 

UNIFORM POWER PROGRAM 


) ! n MI I I j I I[j ! I I M I I I I I I ! I ! ! ! ! j I M ! ! I ! ! I I I I! 11 I M n I I M M I I I I II I I 

2 H ! POWER INPUT PROGRAM ' ' ' 

4 ! ! I M ! I M I ! I I I ! ! ! ! ! ! ! ! ! I I M ! I > ! ! ! I I M I 11 M I ! I ! I M M I I ' ' ' I ' I ' ' I * * ' 

5 !OUTPUT 7091"LOCAL CLR" 

7 OUTPUT 7091*RST" 

10 REAL Uolts<51) 

20 real Poij(5I ) 

30 REAL Ptot ,Pavg .Perdif 

920 OUTPUT 709;"CONFMEAS DCU .308-315.USE 0" 

930 FOR 1-1 TO 8 

931 J-9-I 

940 ENTER 709:Volta(J ) 

9S0 PRINT " “ 

952 NEXT I 

953 OUTPUT 709:"RST" 

954 OUTPUT 709:"CONFMEAS DCU.516.USE 0“ 

955 ENTER 709iUolts(46 ) 

9S£ OUTPUT 709:"RST" 

971 OUTPUT 709:“CONFMEAS DCU .403-417 .USE 0" 

972 FOR I-S TO 25 

973 ENTER 709:Uolt5(I ) 

977 NEXT I 

976 OUTPUT 70S:"RST" 

979 OUTPUT 709: "CONFMEAS DCU .400-401 .USE 0'' 

980 FOR 1-47 TO 48 

961 ENTER 709:Uolts(I ) 

982 NEXT I 

983 OUTPUT 709:“RST" 

985 OUTPUT 709:“CONFMEAS DCU .503-517 .USE 0" 

98B FOR 1-24 TO 36 

987 ENTER 709:Uolt5( I ) 

989 NEXT I 

990 OUTPUT 709;“RST" 

991 OUTPUT 709:“CONFMEAS DCU ,500-501.USE 0" 

992 FOR 1-49 TO 50 

993 ENTER 70g:Uolt5< I ) 

994 NEXT I 

995 OUTPUT 709:"RST“ 

997 OUTPUT 709:"CONFMEAS DCU .502-508 .USE 0" 

998 FOR 1-39 TO 45 

999 ENTER 709:Uolts( I ) 

1000 NEXT I 

1001 OUTPUT 709:"RST" 

1002 OUTPUT 709:“CONFMEAS DCU,600 .USE 0" 

1003 ENTER 709:Uolt5(51) 

1004 OUTPUT 709:"RST” 

1005 PRINT “ " 

1005 PRINT " " 

1050 PRINT "INPUT UOLTAGE D.C. UOLTS PAN 1":Uolts(4G) 

1051 PRINT "INPUT UOLTAGE D.C. UOLTS PAN 2”:Uolts(47 ) 

1052 PRINT "INPUT UOLTAGE D.C. UOLTS PAN 3":Uolts(48) 

1053 PRINT "INPUT UOLTAGE D.C. UOLTS PAN 4":UoHs(49) 


119 






1054 PRINT 'INPUT UOLTAGE D.C. UOLTS PAN 5“;Uolts<50> 

1055 PRINT 'INPUT UOLTAGE D.C. UOLTS PAN 6":Uolt5<5l) 

1060 PRINT “ ■ 

1070 PRINT ■ " 

1080 Resist-2.0 
1108 R2-Resist*Resist 

1119 Pou(2 )-(Uolts(4E )-Uolts(2 ) )*{Uolts(2 )«<Resist + .4 )-Uolts(4B)•.4)/R2 
1149 Pow<5 )-<Uolts(4B )-Uolts(5))«<Uolts<5)*<Resist+.2 )-UoHs(46 )*.2 >/R2 
1179 Poy<9)-(Uolts(47)-Uolts(9))»(Uolts(9>*(Resist+.3)-Uolts(47)*.3)/R2 
1189 Pou( 10)-(Uolts(47 )-Uolts{ 10) )»<UoHs( 10 )*<Reslst+.45 )-Uolt 5 ( 47 )• .45)/R2 
1194 Pou<18)-<Uolts(48 )-Uolts<I8))*<Uolts(18 )*(Resist + .55 )-Uolts<48>•.55 )/R: 

1197 Pow<23)-(Uolts(48)-Uolts<23))«<Uolts<23>*(Resist+,55)-UoIt5(48)*,55)/R2 

1198 Pou<24)-(Uolts(49)-Uolts<24>)»(Uolts(24>*(Re5ist+.1 >-Uolts(48 )*.1 )/R2 

1199 Pow<25 >-<Uolts< 49)-Uolts<25>>*<Uolts(25>*<Resi5t + .3 >-Uolt5<4B >».3 )/R2 

1238 Pow(29)-(Uolts(49)-Uolts(29))»(Uolts(29>«(Resist+.45)-Uolts(49)*.45)/R: 

1241 Pow<32 )-<Uolts(50 )-Uolts(32 ) )*(Uolts<32 )«(Resist+.B)-Uolts(50)*.6)/R2 

1243 Pow<34 )-(Uolts<50)-Uolts<34 > )»(Uolts<34 >«(Resist + .3 )-Uolts(50 )*.3 )/R2 

1244 Poy(35)-(Uolts<50)-Uolts<35)>*(Uolts(35 >*(Re5ist + .B)-Uo:ts<50)».B )/R2 

1245 Pou(3B )-(Uolts(50)-Uolts<36 ))*<Uolts(3B )*(Resist+.6 >-Uolt 5 (50>*.B )/R2 

124G Pou<37)-(Uolts{50)-Uolts<37))*<Uolts(37)«(Re5ist+.B )-Uolts(50 )*.6 )/R2 

1247 Pou<39)-(Uolts<S1 )-Uolts(39))»(Uolt5(39)*(Resi5t + .3)-Uolt5(51 )*.3)/R2 

1248 Pow(40)-(Uolts<51 )-Uolts<40))*(UoltsU0)*(Re5i5t + .3)-Uolts<51 >*.3)/R2 

1250 Pow(42 )-(Uolts<Sl )-Uolt8<42 ) )»(Uolts<42 >»(Resist+.55)-Uolts( 51 )*.55 )/R2 

1251 Pow(43 )-<Uolts(51 )-Uolt5<43>)*(Uolts<43)«(Re5i5t+.45)-Uclts(51 )».45 )/R2 

1252 Pou(44 )-<Uolts(51 ;-UoU 8(44 1 )*(Uolts(44 )*(Resi5t+.55>-Uolt5l51 )*,551/62 
1281 Pou(2 )-(UoIts<46 )-Uolts<2 ) )*(Uolts(2)*(Resist + .4)-Uolts(46)*.4)/R2 
1284 Pou(5)-<UoIt5(4S )-Uolts<5 ))*(yolt5(5>*<Re5i5t+.2 )“Uolts(46 )*.2 )/R2 
1290 Pow<4)-(Uolts<4S)-Uoits<4))*(Uolt5<4>»(Resi5t+.l )-Uolt5(46 )•.1 )/R2 

1294 Pou<9)-<Uolts(47)-Uolt5<9))*(Uolt5(9>*(Re5i5t+.3>-Uolt5(47>*.3)/R2 

1295 Pou(25 )-(Uolts(49)-Uolts(25 ) )*<Uolt5(25 )•(Resi5t+.3)-UoltsC4 9)•.3)/R: 

1297 Pou(29)-<Uolts(4g)-Uolts(29))»(Uolts(29)*(Re5i5t + .2>-Uolt5(49)*.2 )/R2 

1296 Resist-1.972 

1299 R2-Resi5t*Resist 

1300 Poui 15 '-(Uolts(47 .'-Uolt5( 15 ) )*(Uolts( 15 >*(Resist+.2 )-Uolts(47 >«.2 )/R2 

1301 Resist-I.9184 

1302 R2“Re8ist*Resist 

1303 Pou(14 >»(Uolt s(47 )-Uoits<14))»(Uolt5(14)»<Resist+.1 )-Uoltsi47)«.1 >/R2 

1304 Resist-1.886 

1305 R2“Resi5t*Resist 

1306 Pou(4)-<Uolt5(46)-UoIt5<4))*(Uolts(4)*<Rc5i5t+.l )-Uolts(46 )•.1 )/R2 

1307 Pow<6 )-(Uolts(46 )-Uolts<6 ) )•(Uoltsi6 )•<Resxstt.1 )-Uolts<46 )•.1 )/R2 

1308 Pow(30)-(Uolt s(49 )-Uolts(30))*(Uolts<30>»<Resx5t + .4 )-Uolts(49 )*.4 )/R2 

1309 Res X st-1.862 

1310 R2-Resist»Resxst 

1311 Pou(27)-(Uolts(49)-Uolts<27))*<Uolt5(27>«(Resxst+.S)-Uolts(4g)».5)/R2 
'312 Pow(28)-(Uolts(49)-Uolts<28))*(Uolts<28)*(Resxst + 0. )-Uolt 5 (49 )«0. )/R2 

1313 Pou(3l >-(Uolts<50)-Uolts<31 ))«(Uolts(31 )*<Re5xst + .65 )-Uolts(50 )*.65 )/R2 

1314 Pow(38)-(Uolts(50 )-Uolts<38 ) )*(Uolts(38 )•(Resi5t+.3 )-Uolts(50 )•.3 )/R2 
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1315 Resist'l.852 
131S R2“Re5i5t*Resist 

1317 Pow(24)-<Volt5(49)-Uolt3<24 ))«<Uolt5il4)*(Resi5t + ,1)-yolts( 49)*.l )/R2 

1318 Pow(33)-<VoUs(50)-Volts(33))*(Volt5<53>*(Re5i5t+.7)-yolt9(50)*.7>/R2 

1319 Resist*!.8333 

1320 R2*Resist»Resist 

1321 Pou(11 )-(Uolts(47)-yolt5<11 ) )*<Volts<11 )•<Resist*.5l-Woltsi 47)•,5)/R2 

1322 Poij( 12 >-(yolts<47>-Volts( 12 ) )*(yolts( 12 )*<Re5ist+. 1 )-yolt 5 ( 47 )•. ; -'R; 

1323 Pow(I3>»<yolts<47)-yolts(13))*(Volts(13 )*(Resist*.2 l-Wolts(47)*.2)/R2 
1325 Pou(41 )-(Uolts(51 )-Uolts<41 ))*(Uolts<41 )*<Re5ist*,1 )-yolt9(51 )•.) )/R2 
1325 Pow(45 )-(Uolts(51 >-yolts<45 ) )•<Uoltsi45 )*(Resist*.5 l-Oolt5( 51)».6)/R2 

1327 Resist-1.8182 

1328 R2*Resist*Resist 

1329 Pow(21 )-(yolts(48>-yolts(21 ))*<Uolts<21 >*<Resist*. 7 >-yolts(48>♦. 7 )/R2 

1330 Res 1 st»1.8 

1331 R2“Resist»Resi9t 

1332 Pou<8>-(yolts(4G )-Uolts<8 ) )*(Vol19(8 )*(Resist*.2 l-Uol15(46 >*.2)/R2 

1333 Pow(19)-(Uolts(48 )-Uolts<19) )*<yolts(19>•(Resist*.2 l-Wolt 5 (48 >*.2)/R2 

1334 Poij(20>-(Volts<48)-yolts(20))*(yolts(20>*(Re5i5t*.7)-yolts(4S>*,7>/R2 

1335 Pou(22 )-(yolts(48 )-yolt5(22 ) )«(yolts(22 >*(Resist*.2)-yolt 5 < 48)».2)/R2 
1335 Resist*!.6657 

1337 R2*Resist‘Resist 

1338 Pow<1 >-(yolts(45 )-yolts(1 ) >•<yoltsf1 >•(Resist*.3 )-yolts< 45)♦.3 >/R2 

1339 Pow< 3 )-(yolts(45 )-yolts(5 ) )*(yoitsi3 )*(Resist*.7 >-yolts(45 >•.7 )/R2 

1340 Pou( 7 )-< yolt s< 46 )-yol ts(7 ))»<yolt9<7)*(Resist+.2l-yoltsC 45 >*.2)/R2 

134 1 Poy< 16 )■< yolts( 48 )-Volts( 16 ) )*(yolts( 16 )♦(Resist* .4 )-yolts( 48 >• .4)/R2 

1 342 Pou( 17 )-< yolts(48 )-yolt5< 17 ) )*<yolts(17 )*(Re5ist*.4 >-yolis( 48)♦.4 )/R2 

1343 Poy< 26 )■< yolts< 49 l-Volts(25 ) )•(yoJts(26 )*(Resist*.2 >-yoltsi 49)•.2)/R2 

1344 FOR I-l TO 45 

1345 PRINT -BLOCK « "il." POWER • -jPowU).- WATTS- 

1346 WAIT .5 

1347 NEXT I 

1348 Ptot-0 

1349 FOR 1-1 TO 45 

1350 Ptot-Ptot*Poy<I ) 

1351 NEXT I 

1352 Pavg-Ptot/45 

1353 PRINT - ” 

1354 PRINT "AyERAGE POWER - "iPavg,' WATTS" 

1355 PRINT ' ■ 

1356 FOR I-l TO 45 

1357 Perdif-ABS(1Pavg-Poy(I) )/Pavg )«100 

1358 PRINT “HEATER t":!,- PERCENT DEyiATION -;Perdif 

1359 WAIT .5 
136C NEXT I 
1361 END 
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APPENDIX F 

UNIFORM DATA REDUCTION PROGRAM 

10 I * M > M M I ! ! I I ! I M I M ! M ! ' f ! ! > I M ! I M > < I I I I I I I I M I I M I M I I M I I II I I I M I I 

20 I! TEMPERATURE ACQUISITION AND REDUCTION PROGRAM m 

50 M (STEADY STATE) n 

31 !' PROGRAM PERFORMS THE REQUIRED FUNCTIONS NEEDED TO OBTAIN DATA m 

32 !! FROM THE ACQUISITION SYSTEM AND CONVERT THE THERMOCOUPLE OUT- i! 

33 II PUT VOLTAGES INTO TEMPERATURE VALUES. THE TEMPERATURES ARE 'i 

34 !I PRINTED OUT AND THEN REDUCED TO COMPARE WITH THEORETICAL VALUES '' 

40 II I I I I I I I M I I I I I M I I I ! I I I I I I I I I I I I M I I I I ! I I I I I I I I I I I I I I M I II I I I I I I I I I I I 

50 REAL Volt3(G0) 

60 REAL T(60) 

65 REAL Ndt(45).Ndtl(45).Lnu(45),LQr(45) 

66 REAL Pouer,Tanb ,G ,Kpo ,Oelx ,As .Per,Ta.Qconv,Lbar 

67 REAL Of lux,De11 ,Tfi1n ,Beta .Spvol .Dynvi9C .Ocond 

66 REAL K1 nv 1 9C ,K f ,Pr .Theta ,Hx ,Grpi ,Hux 

63 PRINTER IS CRT 

70 PRINT 'ENTER OUTPUT DATA FILE NrME" 

71 INPUT Filenames 

72 PRINT 

73 PRINT "ENTER POUER VALUE IN WATTS" 

74 INPUT Power 

75 PRINT 

76 PRINT "ENTER SHROUD SPACING IN mm" 

77 INPUT Spaces 

76 'PRINT 

79 'CREATE BOAT *THREDi5: .700,1" ,45 

60 'ASSIGN 9Path3 TO "THRED15: .700,1" 

81 'PRINT DATES(TIMEOATE ) 

63 'PRINT Power." WATT INPUT POWER" 

84 'PRINT " 46 STRIPS NO SHROUD" 

85 'PRINT " BOAT FILE THRED15: .700.1" 

66 'PRINT “ “ 

89 I I I I I I I M II I I I I I I I I I I I I I I I I I SURFACE fl-20 

90 OUTPUT 70grCONFMEAS OCV . 100-1 19.USE 0" 

91 FOR I-l TO 20 

92 ENTER 709(Volt9( I ) 

93 T(I)-,0006797 + (25825.1328*Vo1ts(I ))-(607789.2467*(Volts ( I>«Volts( I)))-( 219 

52034.3364. ( Volts(I )"3 ) )+<8370810996.1e74.(Vol15(I>*4))+.039 

94 NEXT I 

95 I M I I I I I I M I I I 1 I I I I I I I . SURFACE 121-140 

96 OUTPUT 709:"CONFMEAS DCV .200-219 .USE 0" 

97 FOR 1-21 TO 40 

90 ENTER 709:Volt5( I ) 

99 T( I)-.0006797 + (25825.1328«Volts(I ) )-<607789.2467*(Volts(I>«Vclts(I)))-(2l9 

52034.3364. (Volt9(I)*3 ) ) + < 8370810996.1874.(Volts<I)'4))+.039 

100 NEXT I 

101 I I I I I I I I I I I ...I SURFAC t41-#45 

102 OUTPUT 709:"CONFMEAS OCV .300-304 .USE 0" 

103 FOR 1-41 TO 45 

104 ENTER 709:Voit9( I ) 

105 T( I >-.0006797.(25825.1328.Volts(I ))-(607789.2467.(Volt8(ILVolts(1)))-( 219 

52034.3364. ( Volt9( I)“3)>+(8370810996.1874.(Volts(I>"4)>+.039 

106 NEXT I 











107 M! PRINT OUTPUT TEMPERATURES 

108 IPRINT " ■ 

109 IFOR 1-1 TO 45 

110 IPRINT "HEATER #"tl,“ TEMPERATURE IN DEG. C -";T(I) 

111 I NEXT I 

112 !OUTPUT BPath3tT(* ) 

115 IPRINT ' " 

116 IPRINT “ BATH TEMPERATURES (TOP TO BOTTOM)" 

117 IPRINT * ■ 

118 OUTPUT 70g!"CONFMEAS DCU.317-31 9 .USE 0" 

119 FOR 1-57 TO 59 

120 ENTER 709iUolts(I ) 

121 T(I)-.0006797+(25825.1328»Uolta{I) >-(607789.2467«<Uolts<I )*Uolt5(I> > )-(219 
52034.3364*<Uolt5(I )*3))+(8370810996.1874*(Uolt5(I )*4) )+.03G 

122 IPRINT "TEMP. DEG. C "|T(I) 

124 NEXT I 

126 T<22)-25 

127 T(29)-25 

133 IPRINT "INPUT AMBIENT TEMPERATURE" 

134 Tanb-(T(57)+T(58)+T(59))/3 

135 PRINTER IS 701 

136 PRINT " " 

137 PRINT DATE$(TIMEDATE) 

158 PRINT "DATA REDUCTION OUTPUT" 

139 PRINT " ” 

140 PRINT "POUER (WATTS) :" .Power 

141 PRINT 

U: PRINT "SHRul- POSITION (mm) Spaces 

143 PRINT 

144 PRINT AMBIENT TEMPERATURE (C) ;" .Tawb 

145 PRINT 

146 PRINT 

147 

148 'ASSIGN 0Fath TO FilenaMeS 

MS 'ENTER 9PathiT(* ) 

150 ' 

151 G-9.81 

1 52 Kpg-.1421 

153 Delx-.006731 

154 A5-.0239*.0078 

155 Per-2«.0239+2*.0078 

156 Lbar-As/Per 

157 A$-"Qcond'' 

158 B$-“X/L" 

159 CS-“NDT" 

160 D$-"L06(Nux )" 

161 E$-"LOG(Grx* )" 

162 FS-"DELT" 

163 GS-"Hx" 

164 H$-"Nux“ 

165 I$-"6rx.*" 

166 J$-"TeMp" 

168 PRINTER IS CRT 
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1E9 I ! ! ! M ! I ! ! ! I ! ! I I ! M I M I > I I ! M ! M ' < I I i M ! I ! ! ! ! ! ! ! 

170 ! ENTER CONDUCTION LOSSES I 

171 ! I ! ! ! H I I I I M M I M I I I M j I I I I ! ! ! I I M ! M M I I ! M M ! 

172 PRINT “INPUT CONDUCTION LOSSES AS COMPUTED FROM ELLPACK” 

173 INPUT Ocond 

174 !I!!! ! ! I ! I ! ! ! ! I ! !! M ! ! !! M I ! ! I ! I M I M I I ! ! > ! I ! !! I 

175 CREATE BOAT FilenaneS .45 

17B ASSIGN #Path3 TO Filenames 

177 PRINTER IS 701 

178 FOR I-l TO 45 

179 PRINT “ " 

180 PRINT • “ 

181 PRINT " SURFACE #“ .1 


182 PRINT * • 

183 PRINT USING “2X.AAAAA,3X.AAA,3X,AAA,3X .AAAAAAAA ,2X .AAAAAAAAA.3X ,AAAA ,4X ,AA 
,4X .AAA .3X .AAAA .2X .AAAA";AS.BS.C$.DS,ES.F$.S$,H$.1$.J$ 

184 Qconv»Pouer-Qcond 

185 Of lux<»Qconv/As 
18G Delt-T(I>-Tamb 

187 Tfilm-273, lS+(Tn >+Tamb)/2 

188 Beta-(-7.9448S8£-8*(Tf ilm'‘2 ) )+<5.735G479E-5»Tf iIn .00978105B3 

189 Spvol-(4.S9S99E-9»<Tfilfn'2 ) )+<-2.5574SE-G»TfiIm )+.001341903 

190 Oynvisc-(3.2548511E-7*(Tfilm“2))+(-.00021474487*TfiIm>+.03G166792 

191 K'invisc“Spvol*Dynvi sc 

192 Kf-(1.4181818l8E-3*Tfiim)+.1866 

193 Pr-(4.65706E-3»<Tfilm*2 ) >+(-2.922094*TfiIm )+463.3319 

194 Thefal’■Power/I .0220*Kf ) 

195 NdtKI )-0eit/Theta1 

196 Nd-Ndtt<I) 

197 Temp-T(I) 

199 Theta*0flux»Lbar/Kf 

200 Ndt( I )-OeU/Theta 

201 Hx“0flux/Delt 

202 6rm-(G*Beta»0flux*<Lbar)*4)/(Kf*Kinvisc“2) 

203 Nux»Hx♦Lbar/Kf 

204 L 9 r(I )«LGT(Grm) 

205 Lnu(I )*LGT(Nux> 

206 Lg-Lgrd) 

207 Lnn«Lnu(I ) 

209 Xtoi-(31.75-(6.985+11-16)*1.270))/31.75 

210 PRINT USING "1X.O.DOOD,2X,D.DDO.IX.D.DDD.2X.D.DDD.5X,D.DDD.4X,DD.DDD,2X,DD 
DD.D .IX .DD.DDD ,DDDDD.D,1X ,DD.00Ocond ,XtoI .Nd .Lnn .Lg .Celt .H.v .Nu.r ,Grm .Temp 

215 NEXT I 

225 OUTPUT 0Path5iT<* ) ,Ndt1(• ) .Ndt(* ) .Lnu(• ) .Lgrl* ) 

226 PRINTER IS CRT 
235 END 
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APPENDIX G 

NONUNIFORM POWER PROGRAM 


2 

4 

!*! POWER 

1 M 1 1 1 M M M M 1 M 1 1 1 

INPUT PROGRAM 

1111 1 u11; 1 ;;;;;; 1 1 111 1 M ; 1 ; 

5 

•OUTPUT 709;“LOCftL CLR“ 

7 

OUTPUT 709rRST" 


10 

REAL UoltsCBl ) 


20 

REAL Pou(51 ) 


30 

REAL Ptot ,PavQ .Perdif 

920 

OUTPUT 709t"CONFMEAS 

DCV .508-315 .USE 0" 

950 

FOR 1-1 TO 8 


951 

J-9-I 


940 

ENTER 709iUolts< J ) 


950 

PRINT " “ 


9EZ 

NEXT I 


955 

OUTPUT 709i'’RST" 


954 

OUTPUT 709i''CONFMEAS 

DCU.316.USE 0“ 

955 

ENTER 709!Uolts(45 ) 


955 

OUTPUT 709;"RST" 


971 

OUTPUT 709;“CONFMEAS 

DCV .403-417.USE 0" 

972 

FOR 1=9 TO 23 


975 

ENTER 709;Uolts(I ) 


977 

NEXT I 


973 

0UT=UT 709;“RST" 


979 

OUTPUT 709;"CONFMEAS 

DCV ,400-401 .USE 0" 

980 

FOR 1-47 TO 48 


951 

ENTER 709;Volts(I ) 


962 

NEXT I 


985 

OUTPUT 709;"RST" 


995 

OUTPUT 709;"CONFMEAS 

DCV .503-517,USE 0" 

986 

FOR 1-24 TO 58 


987 

ENTER 709;Uolt5<I) 


989 

NEXT I 


990 

OUTPUT 709;"RST" 


991 

OUTPUT 709;"CONFMEAS 

DCV ,500-501.USE 0" 

992 

FOR 1-49 TO 50 


993 

ENTER 709;Uolt5{I ) 


994 

NEXT I 


995 

OUTPUT 709;"RST" 


997 

OUTPUT 709;"CONFMEAS 

DCV ,602-508 .USE 0" 

998 

FOR 1-39 TO 45 


999 

ENTER 709;Uolts< I ) 


1000 

NEXT I 


1001 

OUTPUT 709;"RST" 


1002 

OUTPUT 709;"CONFMEAS 

DCV ,600 .USE 0" 

1003 

ENTER 709;Volts<5l) 


1004 

OUTPUT 709;"RST" 


1005 

PRINT •' " 


1005 

PRINT ■ " 


1050 

PRINT "INPUT VOLTAGE 

D.C. VOLTS PAN 1"iVolt5(46 ) 

1051 

PRINT "INPUT VOLTAGE 

D.C. VOLTS PAN 2“;Volts(47) 

1052 

PRINT "INPUT VOLTAGE 

D.C. VOLTS PAN 3";Vol15(48 ) 

10£3 

PRINT "INPUT VOLTAGE 

D.C. VOLTS PAN 4";Volts(49) 
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1054 PRINT 'INPUT UOLTftSE D.C. UOLTS PAN 5"iUolt5(S0) 

1055 PRINT “INPUT UOlTAGE D.C. UOLTS PAN S’lUoltsiSl ) 

1060 PRINT " " 

1070 PRINT " “ 

1080 Resi9t>2.0 
1109 R2*Re9i5t»Re5ist 

1119 Pou(2 )-(Uolts<4S )-Uolts<2 ) )•(Uolts(2 )•(Resistt.4 )-UoltsC46)«.4 >/R2 
1149 Pow(5 )-<Uolts<46 )-Uolts(5 > )♦(Uolts(5)*(Resist +.2)-UoJts(46)*.2 )/R2 
1 179 Pow(9)-(UoH5(47)-Uolts(9) )*(Uolts(9>«<Re9ist+.3>-Uolt9(47 )*.3 )/R2 

1 189 Pou( 10 )-(Uolt5<47)-Uolts( 10 > )• <Uolts( 10 )*(Resist + .45 >-Uolt9<47 ,'».45 )/R2 
I 194 Pow( 18 >-(Uolt5( 48 )-Uolts( 1 8 ) )*< Uol t s< 1 8 )•< Resi 91+. 55 )-Uolt9(48 )-'.55 )/R2 

1197 Pou(23 )«<Uol19(48 )-Uol19(25 ) )*<Uolts(25 )*<Resi91 +.55 >-Uolt9<48 )».55 )/R2 

1198 Poy(24)-(Uolt5(49)-Uolts<24>)*(Uolt9(24)«(Re9X9t+.l )-Uolt9<48>•.1 )/R2 

1199 Pow(25)-(Uolt5(49)-Uolt9(25) >*(Uolt5(25 >•(Resi9t+.3)-Uolts(48 )».3)/R2 

1238 Pow(29)-(Uolts(49)-Uolts(29))«<Uolts(29>«(Re9ist+.45)-Uolt9(49)».45)/R2 

1241 Fow(32 )-<Uolts(50)-Uolts<32))*(Uolts(32)•(Resist+.6)-Uolts(50)*.6)/R2 

1243 Pow(34 )-(Uolts(50 )-Uolts<34 ))•(Uol15(34 )*(Resist+.3 >-Uolts(50 )*.3 )/R2 

1244 Pouj(35 )■( Uol 15(50 )-Uolts(35 ) )»(Uolts(55 )*( Resists .6 >-Uolts( 50 )* .6 )/R2 

1245 Poy(36)»(Uolts(50 )-Uolts(3G ))•(Uol19(36 >«(Re9ist + .6)-Uolts(50)*.6 )/R2 

1246 Pow(37)-(Uolts(50 >-Uo11s(37 ) )•(Uolts(37)*(Resist+.6)-Uo115( 50)•.6)/R2 

1247 Pow(39)-(Uolt5(51 )-Uolts(39 ) )•(Uol15(39 )•(Resist+.3 )-Uol15(51 )».3)/R2 

1248 Pou(40)*(Uolts(51 )-Uolts(40)>*(Uolts(40)»(Re5i5t + .3)-Uolt5(51 )*.3)/R2 

1250 Pou(42 )»(Uolts(51 )-Uolts(42 ) )*(Uolts(42 )«(Re5ist+.55 )-Uolt5(51 )».55 //R2 

1251 Pow(45 )-(Uolt9(51 )-Uolt9(45))*(Uolt5(43)*(Resi5t+.4S)-Uolt5(5l )*.45)/R2 

1252 Pou(44)-(Uolts(51 )-Uolts(44 ) )•(Uo11s(44>*(Resist+.55 )-Uolts(51 )».55)/R2 
1281 Pou(2 )-(Uolts(46 )-Uolts(2 ) )•(Uolts(2)*(Resist+.4)-Uolis(46)•.4)/R2 

1284 Pou(5 >-(Uolts(46 )-Uolt5(5 ) )•(Uolts(S )•(Resist+.2 )-Uo11s(46 )•.2 >/R2 
1290 Pou(4)«(Uolts(4G)-Uolts(4))»(Uolt5(4)«(Re5i5t+.1 )-Uolts(46)».! )/R2 

1294 Pou(9)-(Uolts(47)-Uolt9(9))*(Uolt5(9)*(Resi9t+.3)-Uolt5(47)*.3 )/R2 

1295 Pow(25 )-(Uolt5(49)-Uolts(25))«(Uolts(25)»(Pesist + .3)-Uolts(49)».3)/R2 

1297 Pou(29)-(Uolts(49)-Uolts(2g ) )•(Uolts(29 )•(Resi5t+.2 )-Uol15(49 )•.2 )/R2 

1298 Resi5t-1.972 

1299 R2“Res 1 51‘Res 1 st 

1300 Pou( 15 )•( Uolts( 47 )-Uoits( 15 ) )»(Uolts( 1 5 )•( Resi st + . 2 )-Uol t s( 47 )♦ . 2 )/F.2 

1301 Res 1 st“1.9184 
1502 R2“Res 1 st♦Res 1 st 

1303 Pow(14)«(Uolts(47)-Uolt5( l4))»(Uolt5(14)«(Resi5tt.l )-Uolt5(47)«.l )/R2 

1304 Resi91“1.886 

1305 R2»Res 1 st♦Res 1 st 

1306 Pow(4 )•(Uolt5(46 )-Uclt5(4 ) )*(Uolt5(4 )♦(Resist+.1 )-Uo11s(46 )♦.1 )/R2 

1307 Pou( 6 1-(Uolt5(46 )-Uolts(6 ) )*(Uolt5(G )*(Resi5t + .1 )-Uo]ts(46 )♦.1 )/R2 

1308 Pou(50)"(Uolts(49)-Uolts(30;)^(Uolt8(30)«(Resist+.4)-Uolt5(49)*.4>/R2 

1309 Res 1 st•1.662 

1310 R2“Res 1 91♦Res 1 st 

1311 Pou(27 )-(Uolts(49 )-Uolts(27 ) )♦(Uoltsl27 )♦(Resist♦.5)-Uo11s(49)♦.5)/R2 
1512 Pou( 28)-( Uolt9(49 )-Uolts(28 ) )^(Uolt5(28 )^(Resi5t+0. )-Uolts(49 )/R2 
1315 Pow(31 >■(Uol19(50 )-Uolts(31 ))»(Uolts(31 )♦(Resist+.65 )-Uo11s(50 )♦.65 )/R2 

1314 Pou(38)-(Uolt9(50)-Uolts(38))^(Uolt9(38)^(Resi9t+.5)-Uolt5(50)^.5)/R2 

1315 Resist“1.852 

1316 R2“Resist^Resist 

1317 Pou(24)-(Uolts(49 )-Uolts(24 ))♦(Uolt5(24 >♦(Resists.1 )-Uolts(49)♦.1 )/R2 

1318 Pou(33)-(Uolts<50 )-Uolts<33 ) )•(Uol19(33 )♦(Resist + .7 )-Uolts(50 )^.7 )/R2 

1319 Resist-l.8333 

1 j 2 0 R2“Res 1 st♦Res 1 st 

1321 Pow(n >-(Uolts(47 )-Uolt5(II ) )♦(Uolts(11 )♦(Resist +.5 )-Uolts(47 )♦.5 1/R2 

1322 Pou( 1 2 )-( Uoits( 47 )-Uol ts( 12 .' >• ( Uol15( 1 2 )♦ ( Resis t + . 1 )-Uol ts( 47 )• . 1 )/R2 

1323 Pou( 13)-( Uolts(47)-Uolts(13 ) >♦(Uoltsl13 >♦(Resist+.2)-Uolt5(47)•.2)/R2 

1325 Pow(41 )»(Uolts(51 )-Uolt9(4! ))#(Uolt9(41 )♦(Resists.1 )-Uolt5(51 )♦.! )/R2 

1326 Pou(45)-(Uolts(51 )-Uolt9(45 ) )^(Uolt5(45)^(Re5ist+.6)-Uolts(51 )^.6 >/R2 
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1327 Resist-!.8182 

1328 R2-Re5i5t»Re5i5t 

1329 Pou(2l)-(VoIts(48)-yolt5(2f >)«(Uo2ts(2l )•fResist+ .7 l-folts(48 )•.7 )/R2 

1330 Resist-!.8 

133! R2-Resist«Re5ist 

1332 Pow(8 >-<Volt5<4e )-Volts(8 >)*<Volts(8)*(Resist+.2 )-yolts(46 )*.2 >/R2 

1333 Pow( 19)-< VoltsC 48)-Volts< 19) )»(Uolt5( 19 >*(Re5ist+.2 )-Volt5<48 )* .2 )/R2 

1334 Pow(20)-(Uolts(48)-Volts(20 ) >»(Volts(20 >*<Re5i5t+.7 I-Woltsl48 )*.7 )/R2 

1335 Pow(22 >-(Volts!48 )-Volts< 22))«(Volts(22)«<Re5i5t + .2 l-Volts!48 )♦.2 )/R2 

1336 Resist-!.6687 

1337 R2»Resi5t»Resist 

1338 Pou (! )-(Volts!46 )-Volts! 1))•<Volts!1)•!Resi 5 t+.3 >-Volts!46)•.3 >/R2 

1339 Pow!3 )•!Volts!46 l-Volts!3)>«!Volt5!3)*!Re5i5t+.7 l-Volts!46 )•.7 )/R2 

1340 Pout! 7 )-! Volt5!46 )-Volts!7 ) )*! Volts!7 )*< Resist +. 2 l-Volts! 46 >• .2 )/R2 
134! Pow!!6 )-!Volts!48 )-Volts!16))•!Volts!16 )•!Resist +.4l-Volts! 48)•.4)/R2 

1342 Pow!17 )-!Volts!48 )-Volts!17) )♦!Volts!17 >*!Resist+ .4 >-Volts!481*.4 )/R2 

1343 Pow!26 1-!Volts!49l-Volts!261}«tUolts!26l»!Resist+.2l-Volts! 491*.21/R2 

1344 FOR 1-! TO 45 

1346 PRINT “BLOCK t "il." POWER - “;Pow!Il,‘ WATTS" 

1347 WAIT .5 

1348 NEXT I 

1350 Ptot-0 

1351 FOR I-l TO 15 

1352 Ptot-Ptot+Pcw!I 1 

1353 NEXT I 

1354 FOR I-3I TO 45 

1355 Ptot-Ptot+Pow!I 1 

1356 NEXT I 

1357 Pavgsides=Ptot/30 
1 353 PRINT '■ " 

1359 PRINT "AVERAGE POWER SIDES- "jPavgsides . “ WATTS" 

1360 PRINT “ ■' 

136! Ptot»0 

1362 FOR I-IG TO 30 

1363 Ptot=Ptot+Pow!I 1 

1364 NEXT I 

1365 Pavgcenter“Ptot/15 

1357 PRINT "AVERAGE POWER CENTER- ";Pavgcenter ," WATTS" 

1368 PRINT " " 

1370 PRINT " " 

1371 FOR I-I TO !5 

1372 Perdif-ABS!!Pavgsides-Pou!ID/Pavgsides 1»!00 

1373 PRINT "HEATER »"!l," PERCENT DEVIATION "(Perdif 

1374 WAIT .5 

1375 NEXT I 

1376 FOR 1-16 TO 30 

1377 Perdif-ABS!!Pavgcenter-Pow!I 1l/Pavgcenter!•!00 

1378 PRINT “HEATER «";I," PERCENT DEVIATION ”;Perdif 
1370 WAIT .5 

1380 NEXT I 

138! FOR I-3I TO 45 

1582 Perdif-ABS!!Pavgsides-Pow!ID/Pavgsides 1*100 

1383 PRINT "HEATER «":I,“ PERCENT DEVIATION "iPerdif 

1384 WAIT .5 

1385 NEXT I 

1386 END 
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APPENDIX H 

NONUNIFORM DATA REDUCTION PROGRAM 


10 I ! ! I I ! ! 1 M I ! ! I I M M I ! ! ! ! I !!!I!!!!!! ! ! j ! M!! ! I I M ! ! M ! ! M I M I I I ! I I I MIM 

20 !! TEMPERATURE ACQUISITION AND REDUCTION PR06RAM !! 

30 M (STEADY STATE) !! 

31 M PROGRAM PERFORMS THE REQUIRED FUNCTIONS NEEDED TO OBTAIN DATA !! 

32 !! FROM THE ACQUISITION SYSTEM AND CONUERT THE THERMOCOUPLE OUT- I! 

33 !! PUT VOLTAGES INTO TEMPERATURE VALUES. THE TEMPERATURES ARE !! 

34 !! PRINTED OUT AND THEN REDUCED TO COMPARE WITH THEORETICAL VALUES I' 

40 I ! I I I ! ! j ! ji I M ! I ! M ! ! ! ! ! ! < M I I > ! ! ! I M ! M ! ! I I !!!!M I!I ! IM * I I I I I I I I I I IM 

50 REAL Volts(G0) 

G0 REAL T<B0) 

B5 REAL Ndt(45) .Ndt1(45>,Lnu(45).L0P(45> 

BE REAL Powerl,Pouer2 ,Tamb ,G ,Kpg ,OeIx .As,Per,Ta.Oconv,Lbar 

B7 REAL Of lux.Delt ,Tf 1 Im .Bet a .Spvol.Dynvisc.Qcond 
68 REAL Kinvisc .Kf ,Pr .Theta .Hx .Grw.Hux 
B9 PRINTER IS CRT 

70 PRINT -ENTER OUTPUT DATA FILE NAME” 

71 INPUT Filenames 

72 PRINT 

75 PRINT -ENTER TWO POWER VALUE IN WATTS” 

74 INPUT Powerl .PowerZ 

75 PRINT 

76 PRINT -ENTER SHROUD SPACING IN mm” 

77 INPUT Spaces 

78 1PRINT 

79 ICREATE BOAT "THREOI5: .700.1”.45 

80 'ASSIGN ePathS TO “THRED15: .700.1" 

81 IPRINT OATESITIMEDATE ) 

83 'PRINT Power.” WATT INPUT POWER" 

84 IPRINT - 45 STRIPS NO SHROUD- 

85 IPRINT • BOAT FILE THREO15: .700.1" 

86 IPRINT - • 

89 I I I I! ! I M ! I I ! I I I I I I I I I I I I I I I SURFACE 11-20 

90 OUTPUT 709;-CONFMEAS DCV . 100-119.USE 0” 

91 FOR I»1 TO 20 

92 ENTER 70g;Volts(I ) 

93 T(I)-.000G797+(25825.1328*Volts(I ) )-<B07789.2467«<Volt5(I)*VoIts(I)))-(2IS 
52034.53G4*(Volts(I)”3 )) + (85708)0996.1874*(Volts(I>”4)>+.039 

94 NEXT I 

95 I I M I I ( ! I I! I I I I I I I ! ! ! I ! j 1 I I SURFACE f21-140 

96 OUTPUT 709;-CONFMEAS DCV.200-219.USE 0” 

97 FOR 1-21 TO 40 

98 ENTER 709;Volts<I ) 

99 T< I)-.000G797+< 25825.1328*Uolt5(I) >-<607789.24B7*(Volls( I)«UoltsiI>)>-( 219 

52034.3364»(Volt s' I>”3 >) + (8370810996.1874*(VoltsCI )“4) )+.039 

100 NEXT I 

101 I!I< M HI I I I I M ! I I I I I I I I I > SURFAC *41-S45 


102 OUTPUT 709;"CONFMEAS DCV .300-304.USE 0" 

103 FOR 1-41 TO 45 

104 ENTER 70g;Volts<I) 

105 T( I)-.000G797+(25825.1328»Volts<I) )-<B07789.2467*(Volts(I>*Voit5(I)>>-(219 
52034.3364*<Volts(I>'3 > >+<8370810996.1874*(Volt si I>'4)>+.039 

106 NEXT I 
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107 !!! PRINT OUTPUT TEMPERATURES 

108 'PRINT • " 

109 IFOR I-I TO 45 

110 'PRINT -HEATER fjl.- TEMPERATURE IN DEG. C -“;T(I) 

111 INEXT I 

112 'OUTPUT 9Path3«T(* ) 

115 IPRINT ■ - 

116 IPRINT • BATH TEMPERATURES (TOP TO BOTTOM )" 

117 'PRINT • “ 

118 OUTPUT 709|-CONFMEAS DCU .317-319 .USE 0" 

119 FOR 1-57 TO 59 

120 ENTER 709;Uolt3<I) 

121 T<I>-.00067g7+(25825.1328*Uqlts(I) )-<607789.2467«(Uolts(I )*Uolts(I )) )-<2l9 
52034.3364»(Uolt5( I )''3 ) )+<837081099B.. 1874*(Uolts( I >“4 ) )+.036 

122 IPRINT -TEMP. DEG. C “tT<I> 

124 NEXT I 

125 T<22)-25 

126 T(29)-25 

133 IPRINT “INPUT AMBIENT TEMPERATURE- 

134 Tanb-(T(57)+T<58>+T(59> )/3 

155 PRINTER IS 701 

136 PRINT “ “ 

157 PRINT DATES'TIMEDATE ) 

138 PRINT "DATA REDUCTION OUTPUT- 

139 PRINT ” " 

140 PRINT -POWER LEUELS (WATTS) .Power 1 .PouerZ 

141 PRINT 

142 PRINT -SHROUD POSITION (nn) .Spaces 

U3 PRINT 

144 PRINT -AMBIENT TEMPERATURE (C) :”.Tanb 

145 PRINT 

146 PRINT 

147 ' 

148 'ASSIGN @Path TO Filenames 

149 lENTER 9PathiT(* )• 

150 I 

IS! S-9.81 

152 Kpg-.1421 

155 Oelx-.006731 

154 As-.0239*.0078 

155 Per=2*.0259+2».0078 

156 Lbar-As/Per 

157 A5="0cond" 

159 B$-"X/L" 

159 CS--NDT" 

160 DS--UOG(Nux )“ 

161 E$-“LOG(Srx* )- 

162 F$--DELT- 

163 GS--Hx“ 

164 HS-”Nux" 

165 I$-“Grx*’ 

166 JS-'Tenp- 

168 PRINTER IS CRT 

153 I I I I M I I n M M I I I u M 11 I M I I M M ' ' ' I I I I ' I I ' I I > ' I 

170 I ENTER CONDUCTION LOSSES I 

171 M I M I M I I I I I I ! I M M I M I M I I M I M I M I I I I M n I I M 

172 PRINT “INPUT SIDE COLUMN CONDUCTION LOSSES AS COMPUTED FROM ELLPACK- 

173 INPUT Ocond 

174 I I I I M I n I I I II M M I M I I I M M I ! ! I I M M I I I I M I I n I 

175 CREATE BOAT FilenameS ,45 

176 ASSIGN •Path3 TO Filenames 

177 PRINTER IS 701 
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SURFACE f.I 


«78 FOR >1 TO 15 

179 PRINT • • 

180 PRINT " * 

181 PRINT ‘ 

182 PRINT • ” 

183 PRINT USIN6 "2X ,AAfWPf"i5f .AAA .3X ,AAA ,3X .AAAAAAAA ,2X .AAAAAAAAA ,3X .AAAA ,4X .AA 
.4X .AAA .3X .AAAA.2X.AAAA”iA$.8$.C$.0*.E$.F$.GS.H$.1$.JS 

184 Qconv“Power1-0cond 
IBS Oflux“0conv/A5 

I8G Delt«T<I )-Tamb 

187 TfilM-273.15+(T(I)+Tamb)/2 

188 Beta-(-7.9448S8E-8*(Tfiim"2))+(5.735B479E-5*Tfiln>-,0097810563 

189 Spvol-(4.69B99E-9*<Tfilm*2>>+<-2.S3745E-6*Tfiln)+.00l341903 

190 Dynvi5c-(3.2348511E-7*(Tf»ln“2 ))+(-.0002U74487*Tfil«)+.0361B6792 

191 Kinviac«Spvol*Dynvisc 

192 Kf-(1.418181818E-3*Tfiln )+.18S6 

193 Pr-(4.6S70BE-3*(Tfiln''2))+<-2.922094*Tfiln)+463.3319 

194 Theta1-Pouert/<.0220»Kf> 

195 NdtKI )-Delt/Thetal 

196 Nd-Ndtl(I) 

197 Temp=T<I) 

198 Theta»Qf lux*Lbar/K'f 

199 Ndt(I)-Delt/Theta 

200 Hx-Qflux/Delt 

201 6rfn»( 6«Beta»Qf lux» ( Lbar)''4)/<Kf'»Kinvi5c''2> 

202 Nux«Hx«Lbar/Kf 

205 Lgr<I )”LGT(Grn) 

204 Lnu(I )”L6T<Nux) 

205 Lnn»Lnu(1) 

206 Lo»Lgr(I) 

208 Xtol-(51.75-<G.98S+(1-16)•1.270))/3! .75 

209 PRINT USING “1X ,D.0000 .2X .0.000 , tX.0.ODD .2X ,D.ODD.5X.D.ODD,4X.DD.ODD.2X,DD 
OO.D,1X,00.000,00000.0,1X ,00.00"iQcond,Xtol ,Nd ,Lnn ,Lq .Delt ,Hx,Nux,Grm ,Tenp 

210 NEXT I 

211 PRINTER IS CRT 

217 > M ! I ! ! M > ' M I M M I I! I I ! I I M I ! I I I I! If ! I I I I I M M I 

216 I ENTER CONDUCTION LOSSES ! 

219 ! M M ! I M 1 M M M I If ! M I I ! ! ! I < ' M I I ! I I ! ! I! M ! I M I 

220 PRINT "INPUT CENTER COLUMN CONDUCTION LOSSES AS COMPUTED FROM ELLFACK" 

221 INPUT Ocond 

222 ! M I M I M M M M I M I I U I I I I n I I I I n I M M I I M I I I I I 

225 PRINTER IS 701 

226 FOR 1-16 TO 30 

227 PRINT " “ 

228 PRINT " " 

229 PRINT " SURFACE #".1 

230 PRINT " ■ 

231 PRINT USING "2X .AAAAA .3X.AAA ,3X,AAA.3X.AAAAAAAA .:x.AAAnAAAAA.3X.AAAA ,4X,AA 
,4X,AAA .3X .AAAA ,2X .AAAA"|A$ .B$ ,C* .05 .E$,FS,GS .H$ ,1$,JS 

232 Qconv-Power2-0cond 

233 Oflux-Oconv/As 

234 Delt-T(I )-Tanb 

235 Tfilrt-273.15+<T(I )+Tanb )/2 

236 Beta-(-7.g44858E-8*<TfiIn''2 ))+(5.735B479E-5*Tf iln >-.0097810563 

237 Spvol-<4.B9699E-g*(Tfiln"2) »+<-2.S3745E-6*Tfiln)+.00134)g03 

238 Dynvj5c-(3.234851lE-7*<Tfiln"2))+(-.00021474487*Tfiln>+.036166792 

239 Kinvisc-Spvol*Oynvisc 

240 Kf-(1.418l81818E-3*Tfiln >+.1866 

241 Pr-<4.B570BE-3*(Tfiln'2))+<-2.922094»TfiIn 1+463.3319 

242 Thetal-Power2/(.0220»Kf) 
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243 Ndt1<I)-Delt/Theta) 

244 Nd-Ndtt<l) 

245 Tenp«T <I> 

246 Theta-Qflux*Lbar/Kf 

247 Ndt<I)-Delt/Theta 

248 Hx»Oflux/Delt 

249 G'-n-<e*Beta*Qflux*(Lbar)*4)/(Kf*Kinvjac"2) 

250 Nux-Hx*Lbar/Kf 

252 LgrCD-LGKGrM) 

253 Lnu(I )«LGT(Nux ) 

254 Lnn»Lnu(I) 

255 Lg“LQr<I) 

255 Xtol-(31 .75-(5.98S+< I-16 )•■! .270) )/3l .75 

257 PRINT USING " 1X.0.DODD.2X.0.DOO.1X.D.DDD.2X.0.ODD,5X,D.ODD,4X.DO.ODD ,2X .DD 
DD.O.IX .DD.DDD.DD0DD.D.1X,OD.OD";Ocond.Xtol.Nd.Lnn.Lg.Delt,Hx,Nux.Grn.Tenp 

258 NEXT I 

259 PRINTER IS CRT 

260 ! ! ! I M !! ! I M ! ! ! ! ! ! M ! I ! I ! M I ! I < ! M ! ! ! ! ! !! !!H > ! ! 

261 ! ENTER CONDUCTION LOSSES ! 

262 ' I ! M M I !! ! I! ! I I j I > I I I ! ! ! ! ! M I I M ! ! ! ij ! ! ! ! ( I M ! I 

263 PRINT “INPUT SIDE COLUMN CONDUCTION LOSSES AS COMPUTED FROM ELLPACK' 

264 INPUT Qcond 

265 ! I ! ! M ! M ! i ! ! ! ! ! M ! ! ! ! ! ! I >1 I ! I I ! I I M ! ! I ! I ! ! M ! ! ! 

266 PRINTER IS 701 

267 FOR 1-31 TO 45 

268 PRINT " * 

269 PRINT “ " 

270 PRINT " SURFACE t“.1 

271 PRINT " “ 

272 PRINT USING “2X.AAAAA.3X.AAA,3X.AAA,3X.AAAAAAAA,2X.AAAAAAAAA,3X .AAAA ,4X ,AA 
.4X,AAA ,3X ,AAAA,2X.AAAA"iA$.8$,C$.0$,ES.F$.6$,HS,1$.J$ 

273 Qconv-Pouerl-Ocond 

274 Oflux-Oconv/As 

275 Delt-T(I)-Tamb 

276 Tf jIm- 275.15+<T(I)+Tanb)/2 

277 Beta-(-7.944858E-8*(Tf ilM''2 ))+(S.735G479E-5*TfiIn >-,0097810563 

27S Spvol-i4.69699E-9*(Tfiln“2))+(-2.53745E-6*TfiItn >+.001341903 

279 Dynvi5c-(3.2348511E-7»(Tfiln'2 ))+(-.00021474487*TfiIm >+.036166792 
2S0 Kinvi9c=5pvol»0ynvi5C < 

261 Kf-(1.418181818E-3*Tfiln)+.1866 

282 Pr-(4.6570GE-3*(Tfiln“2 > >+<-2.922094»Tfiin>+463.3319 
263 Theta 1“Power 1/(.022O»Kf ) 

284 Ndt1 I I>-Deli/Thetal 

285 Nd-Ndtl(I> 

286 Tenp-T(I> 

287 Theta-Qflux»Lbar/Kf 

288 Ndt(I )-Delt/Theta 

289 Hx-Oflux/Delt 

290 Grn-( S*Beta»0flux*(Lbar>"4>/(Kf»Kinvi5c''2> 

291 NuA-Hx*Lbar/Kf 

292 Lgr(I>-L&T(Grn> 

293 Lnud >-L6T(Nux > 

294 Lnn-LnuiI> 

295 Lg“Lgr<I> 

297 Xtol-(31.75-(6.98S+(1-16 >*1.270> >/31.75 

298 PRINT USING “IX.D.DDDD,2X,D.DOD,IX.D.DDD,2X.D.DDD,5X,D.DDD.4X.DD.ODD .2X.DD 
DD.D,IX.DD.DDD.DDOOD.D.IX.DD.DD"lOcond.Xtol,Nd .Lnn.Lg.Delt ,Hx.Nux.Grn .Tenp 

259 NEXT I 

300 OUTPUT •Path3;T(*>.Ndt1(*>.Ndt<*>.Lnu(»),Lgr<«> 

301 PRINTER IS CRT 

302 END 


\ 

( 


i 
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